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Abstract 



Existence, uniqueness and stability of solutions is studied for a set of nonlin- 
ear fixed point equations which define self-consistent hydrostatic equilibria of 
a classical continuum fluid that is confined inside a container A C and in 
contact with either a heat and a matter reservoir, or just a heat reservoir. The 
local thermodynamics is furnished by the statistical mechanics of a system 
of hard balls, in the approximation of Carnahan-Starling. The fluid's local 
chemical potential per particle at r G A is the sum of the matter reservoir's 
contribution and a self contribution —{V* p)(r), where p is the fluid den- 
sity function and V a non-negative linear combination of the Newton kernel 
VN(|r|) = — |r|~^, the Yukawa kernel Vydrl) = — |r|~-'-e~'*l""l, and a van der 
Waals kernel l^,v(|r|) = —(1 + x^|r|^)^'^. The fixed point equations involving 
the Yukawa and Newton kernels arc equivalent to semilinear elliptic PDEs of 
second order with a nonlinear, nonlocal boundary condition. We prove the 
existence of a grand canonical phase transition, and of a petit canonical phase 
transition which is embedded in the former. The proofs suggest that, except 
for boundary layers, the grand canonical transition is of the type "all gas -H- 
all liquid" while the petit canonical one is of the type "all vapor -H- liquid drop 
with vapor atmosphere." The latter proof in particular suggests the existence 
of solutions with interface structure which compromise between the all-liquid 
and all-gas density solutions. 
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1 INTRODUCTION 



The interface between physically coexisting thermodynamic (locally) pure phases 
poses a challenging array of problems in statistical mechanics which fall somewhere 
inbetween the micro- and macroscopic realms. The large scale (macroscopic) geome- 
try of the interface can be successfully modeled as a sharp Gibbs interface, computed 
from some constrained principle of minimum surface area (Wulf shape); this gener- 
alizes to the dynamical domain in form of motion by mean curvature and related 
principles. The transversal structure of the interface, which obviously is not resolved 
when the interface is modeled as a Gibbs interface, is the hard problem that lives at 
the fringe of the macro- world and for which there is no definitive answer yet. 

To get a hand on the transversal structure it is customary to invoke a van der 
Waals (for fluids) or Weiss (for magnets) mean-field approximation which allows 
one to study both the large scale geometry and the transversal interface struc- 
ture. While this approximation overly simplifies the problem, it is far from be- 
ing understood completely and continues to attract the attention of mathematical 
physicists. ^ ^ ^ 1^ ^ For the liquid- vapor interface at equilibrium, which is 
the motivation for this paper, the van der Waals type models emerge in Kacl^ type 
scaling limits from the statistical mechanics of systems of interacting microscopic 
particles, with particle densities resolved on the long distance scale of the attractive 
part, Va, of the particle interaction Vr + Va, while the short distance repulsive part 
Vr is absorbed into the local thermodynamics. ^ 

The local thermodynamics of an iV-body system with repulsive pair interaction 
Vr is given by a pure phase, defined in the thermodynamic limit of a macroscop- 
ically spatially uniform system in thermal equilibrium with a heat reservoir at re- 
ciprocal temperature f3 G and a matter reservoir at logarithmic fugacity (i.e., 
chemical potential per particle : temperature ratio) 7 G M, characterized by a 
position-independent pressure : temperature ratio p = p{f3,'-y) and particle density 
f] = d^p{(3, 7) at all points of differentiability of 7 p{(3, 7). On general thermody- 
namic grounds, (/3,7) ^ p{l3,'y) is strictly positive, increasing in 7 G M, and convex 
in /3 G M.^ and 7. By convexity, 7 p(/3,7) is differentiable a.e., but the models 
from physics are expected to be better behaved and feature only finitely many points 
of non-differentiability, at 71 (/3), 72(/3), say. At such a 7fc(/3) typically two differ- 
ent pure phases are equally likely, one of them denser than the other, and one needs 
to select the one which furnishes the local thermodynamics. 

In this paper the local thermodynamics is chosen to represent a continuum formed 
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by many identical hard microscopic balls, known (in a fluid state) as a hard-sphere 
fluid and more generally as a hard-sphere system. A hard-sphere system is charac- 
terized by a /3-independent pressure : temperature ratio, i.e. p{(3,j) = p,{"f). We 
will write ^^(7) for d^p,{'~f). The function ^,(7) has a point of non-differentiability 
at 7fe associated with a fluid- versus-solid transition. Here we are interested in study- 
ing the fluid phases, but for our investigations we do need to have control over this 
singularity. 

Physically, a hard-sphere fluid may model the short distance repulsion between 
the spherical atoms in a noble gas or between neutrons in a neutron fluid. Over 
somewhat larger distances r any two such atoms or neutrons also feel attractive 
forces —Wa, the van der Waals (Jr.) force in the case of atoms, which is due to self- 
induced dipole-dipole interactions associated with their flrst excited conflgurations, 
and the Yukawa force in the case of neutrons, which is explained in terms of the 
pion exchange of the strong nuclear forces. When the number of atoms or neutrons 
becomes too large, as in (helium) brown dwarf stars or in neutron stars, Newtonian 
gravity has to be added. We choose the Va interaction to mimic any of these physical 
systems. More precisely, writing aV for Va, we let aV stand for any non-negative- 
linear combination of the form 



aV{r) = AwKv(r) + A^V^{r) + A^V^r), (1) 

where 

V^{r)= -(l + xV)-^ (2) 

V^{r)= -e-'^Vr, (3) 

V;(r)= -1/r, (4) 

are integral kernels of strictly negative definite compact operators on L^{A) for any 
bounded A C M^, and where G {O.a^v}. Ay G {0,aY}, and Aj^ G {O,^^,}, while 
ctw, C(y, and are strictly positive coupUng constant : temperature ratios. In the van 
der Waals approximation the effect of aV on the system is accounted for by adding 
to the externally generated chemical potential per particle : temperature ratio 7 the 
chemical self potential per particle : temperature ratio at r, given by —{aV * r])^{r), 
where 

{V*7j\{r)= [ V{\r-rMr)d'f. (5) 
Ja 
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We refer to (Vj^ * '7)a(i") the Newton -, to (K- * ^)a('^) Yukawa and to 

(Vw * ^)A(r) as the van der Waals potential of rj at r. 

In liis original study, van der Waals assumed boundary effects to be negligible 
and the density function rjir) to be spatially uniform, i.e. ?7(r) = rj. These assump- 
tions can be rigorously correct only in the infinite volume limit when the fluid 
fills all space uniformly, with gravity "switched off;" note that Vw(| ■ |) ^ L^(R^) 
and V^{\ ■ I) G L1(M3), while ■ |) e L\^M^) merely. When = Q and the 
constant function ?7(r) = r/ is substituted in ([5]) with A = R^, then —{V* rj)^^ = 
v\\Vi\ ■ l)ll^i(.3) is a constant function, too. Setting \\V{\ ■ \)\\^^^^^^ =■ \\V\\^ for 
short, the van der Waals densities fj^^^ are then computed from the van der Waals 
fixed point equatior]0 

V = p'.ii + (^\\y\\,v)- (6) 

In (a, 7)-parameter space there are disjoint, open two-dimensional domains where 
the algebraic equation has one or three solutions in the fluid density regime, 
respectively (see also sections IV & V); these regions are separated by a closed one- 
dimensional subset featuring two solutions of except for one point (the critical 
point) at which only one solution exists. Constant (large enough) a sections and 
constant (intermediate size) 7 sections through the fluid solution manifold over the 
{a, 7) half plane each produce an ^-shaped curve associated with the famous "van 
der Waals loop." In the region with three fluid density solutions, the largest solu- 
tion is interpreted as the liquid density phase, the smallest as the gas (a.k.a. vapor) 
density phase of the fluid, and the intermediate density solution as a thermodynam- 
ically unstable artifact of the model. The liquid and the gas density solutions are 
each stable fixed points of (jS]) under iteration, the intermediate density solution is 
not. However, thermodynamically the liquid and gas density solutions are simulta- 
neously stable only along the gas & liquid coexistence curve a (-)■ 7 = 7g^*(a) of the 
model, determined by Maxwell's equal-areas construction, rigorously vindicated 
in Ref . [25] , while away from this curve (still in the three-solutions region) thermody- 
namically only one of these two solutions is stable, the other one at most metastable. 
Here, thermodynamic stability and metastability are understood with a and 7 fixed, 
and explained below. 

-^In textbooks (e.g. [23], [38]) one usually finds discussions of (|6]) with ^^(7) replaced by van der 
Waals' ptdw(7) which corresponds to a system of many hard rods on a line. While for systems of 
hard balls it gives quantitatively wrong answers, qualitatively they reproduce those obtained with 
the correct p',{j)- Also, usually a value for ||V||li is given without specifying V. 
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More interesting than (E]) is the nonhnear fixed point problem 



Vir) = p'.(7 - (aV * vX^r)) (7) 

in the positive cone of the non-separable Banach space C°(M'^) of bounded con- 
tinuous functions ^^(r), r G M^. If V G L^(M^), then ([7]) can be solved with the 
Ansatz ri{r) = r/, which leads back to the algebraic fixed point equation ([6]). Yet, 
for a hard-sphere fiuid with ^(|r|) = V^drl), whenever (a, 7) is a point on the 
gas & liquid coexistence curve of locally uniform phases computed from ([HD, then 
(modulo translations and rotations) a unique monotonic planar interface solution 
ri{r) = rj{x) exists, where x G M is a cartesian coordinate of M^. This can be 
shown by adapting the ODE arguments on p. 40-41 of Ref . [42] . which are avail- 
able because (—A + K^)VY(|r|) = — 47r5(r). A monotonic planar interface solution 
illustrates the physical phenomenon of coexistence of the liquid and the gas density 
phases; they have been extensively studied in one dimensional models. ^ ^ 

Furthermore, using the equivalent radial ODE problem obtained with the help of 
(—A + K^)VY(|r|) = — 47r(5(r), Mironescut^S] j^as shown that solutions in with 
spherical droplet / bubble geometry exist; these solutions do not exist exactly on the 
gas & liquid coexistence curve for the uniform phases, yet are nearby. Such ODE 
arguments are not available for a hard-sphere fiuid with l^(|r|) = Vw(|r|), and the 
existence and classification of the non-constant solutions in M.^ of ([7]) in this case is 
largely unexplored territory. We also note that since ■ |) ^ L^(]R^), the fixed point 
problem ([7]) is not well defined in as it stands with * t] given by ([5]); however, 
replacing * 77 by 0n and stipulating the familiar Poisson equation A0n = 4:717], 
solutions in for the related PDE problem A0n = 47rp',(7 — a0) do exist; it is 
easy to numerically compute radial solutions, which have applications in planetary 
science, 122] ^o summarize, non-uniform van der Waals fiuid theory furnishes an 
accessible model to study the structure of non- uniform density functions //(r) of a 
hard-sphere fiuid in M.^. Evidently, boundary effects are absent in M^. Moreover, for 
V = Vy and V = V,^ simple ODE techniques greatly facilitate the computation of 
solutions in M^. 

Beside the structure of interfaces, their fiuctuations are of interest. Unfortunately, 
in unbounded space all interface solutions are thermodynamically neither stable 
nor metastable and interface fiuctuations diverge, ^^^ ^"^ ^ even though droplets 
may be quite long-lived in dynamical calculations with the Alan-Cahn and related 
evolution equations. To obtain finite fiuctuations one needs to stabilize the interfaces. 
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An intuitively obvious way to obtain a stable fluid interface is to enclose the 
fluid inside a container A (with either wetting, non-wetting, or neutral mechanical 
boundary conditions) and to replace the thermodynamic contact condition of pre- 
scribed logarithmic fugacity 7 by the stricter one of prescribed amount of matter 
J^r]{r)d^r = N. When A is macroscopic and N halfway inbetween the values of |A|r/ 
for the large and small fluid density values fj solving (E]), then there is too much 
matter in the container to be all vapor, and too little to be all liquid. In this case 
the system must find a compromise structure: either a drop of liquid surrounded 
by vapor or a bubble of vapor inside liquid, depending on the mechanical boundary 
conditions. It is reassuring to find this scenario confirmed numerically for V = Vyf 
and neutral mechanical boundary conditions, ^^31 [32] gj^^ particle simulations of 
many hard balls with attractive — interactions. Interestingly enough, in 

Ref . it was found numerically that the thermodynamic transition from the vapor 
state to the liquid-drop state is not gradual but occurs at a petit-canonical first-order 
phase transition which is embedded in the grand-canonical first order phase transi- 
tion between vapor and liquid. To rigorously prove this empirical picture correct is 
an interesting mathematical problem which is still largely open. 

To make a modest contribution toward its solution we here continue our previous 
[00] ^ 

studyl^ where we numerically evaluated the fixed point problem 



for a cohesive "hard-sphere continuum" inside a container A C with neutral 
boundary. In this paper we study the existence, uniqueness, structure, and stability 
of fluid solutions to (|H]) with rigorous analysis. Stability is defined as follows. 

The stability of solutions of ([H]) for the thermodynamic contact conditions "heat 
and matter reservoirs" is determined by the functional 



which we rigorously derived from the grand-canonical ensemble in Ref. [33] . Solutions 
of ([8]) are critical points of (|9]) in the positive cone of the separable Banach space 
C°(A). A solution r]^ of is globally V stable if "P^^^K] = PA(^^l)^ where 





nnVl\= [ p.il-iaV*r]\ir))d'r + l [ [ 



aV{\r — r\)r](r)r]{v) d^r d^r, 



(9) 





(10) 
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global maximizers are denotecQ ?7^'^(r), their dependence on a, 7 implied. A solution 
Va of & is locally V stable if 

^a^'>,^)l,, <0 (11) 

for all a ^ such that < 77^ + (x G C^{A). Here, 



\ I P'iil - (cyV * r^J^(r)) (aV * a)l (r) d'r 

+ - I I aV{\r -t\)a{r)(r{i)d^rd^f (12) 
2 ./a ./a 



is the diagonal part of the second Gateaux derivative of at r^^. A solution t]^ of 
dH]) satisfying ([n]) but not is called P metastable. If P^/^(a, a)|^^ > ( = ) for 
at least one a, then 77^ is called P unstable {locally V indifferent) . 

In Ref . |33] we also explained that stability of solutions of (|8]) for a given amount 
of matter in thermodynamic contact with a "heat reservoir" at inverse temperature 
/3 (oc a) is defined in terms of the following thermodynamic free energy functional. 
For each density function i] G C°(A) we define: 

(i) its amount of matter in A, 

Af^[q] = I n{r)d^r; (13) 
Ja 

(ii) its energy : temperature ratio, 

^alv] = l^^iv] + 1 f f "^(k - f|)r?(r)r/(f) d^rd^f; (14) 

J A J A 

(iii) its (strictly) classical entropy, 

r]{r)\nr]{r)d^r - r]{r) / {p,{rf) -rf) /rf dlq d^r, (15) 
Ja Jo 

where p,ijf) is the hard-sphere pressure : temperature ratio as function of r^; 

(iv) its free energy : temperature ratio, 

^M=^M'S:;[r,]. (16) 



^"GC" stands for grand canonical because those densities comprise the support of the grand 
canonical measure in the van der Waals limit. 
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The thermodynamic free energy : temperature ratio F^{a,N) is then given by 

F^ia,N) = mm{J^^[ri] \ M^[ri] = N}. (17) 

V 

Solutions of dH]) which saturate ( IT7|) are called globally J-' stable and denotec^l r7^'^(r), 
their dependence on a, N implied. Local T stability etc. is defined in terms of the 
diagonal part of the second Gateaux derivative of J-"^, 

= -l /"<(^A(r))^'(r)rfV + ^ [ [ aV{\r-i\)a{r)a{i)d^rd^f (18) 

where 5,(77) is defined by S^[ri] = fj^s,{ri{r))d^r. Variation is to be carried out 
under the constraint ad^r = so that the r^^-disturbances preserve the number of 
particles. 

We close this introduction by re-emphasizing that our study of the finite volume 
fixed point problem (|8]) is not motivated by trying to understand so-called finite- 
size corrections to dominant infinite volume results. Rather it is motivated by the 
physical phenomenon of stable interface solutions in finite containers holding a fixed 
amount of fluid. What makes such a study difficult are the following two points: (i) 
boundary layer effects are as big as interface effects, and one has to separate these 
in the analysis; (ii) one has to rule out competing solutions which take values in the 
crystal regime, and this leaves little wiggle room in parameter space. To level the 
ground we first study the simpler problem when the amount of fluid is controlled by 
a matter reservoir before turning to the problem with a fixed amount of fluid. 

The rest of this article is structured as follows: 

• In section II we define the Carnahan-Starling approximation to the fluid part 
and the Speedy approximation to the solid part of the function 7 1— )■ ^,(7). 

• In section III we identify a parameter region in the {a, 7) (half) space in which 
all solutions of (|H]) take values exclusively in the fluid density range. 

• In section IV we identify a region in which fluid solutions are unique. 

• In section V we identify a region where various fluid solutions exists. 

• In section VI we study the thermodynamic stability of the fluid solutions in 
contact with heat-plus-matter reservoirs. We prove the existence of a vapor 4-> liquid 
phase transition in the finite- volume grand-canonical ensemble. 

3iipQ)! stands for petit canonical because those densities comprise the support of the petit 
canonical measure in the van der Waals limit. 
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• In section VII we address the thermodynamic stability of the fluid solutions in 
contact only with a heat reservoir and also explain the relationship with the infinite- 
volume Lebowitz-Penrose results. We prove the existence of a petit-canonical 
finite-volume phase transition between the vapor state and a state for which we 
present evidence that it is of liquid-drop type. 

• Appendix A supplies some explicit evaluations valid for spherical geometry. 

• Appendix B lists the nonlinear partial differential equations associated with our 
fixed point integral equations. 

• Appendix C provides a "dictionary" to translate our dimensionless notation 
into conventional dimensional notation used in the physics literature. 

• Appendix D contains a brief list of (minor) errata for our previous papers on 
the subject, Ref.[29] and Ref . |33] . 



2 THERMODYNAMICS OF HARD SPHERE SYSTEMS 

Numerical simulations ^2S] of the dynamics of many identical hard balls indicate 
a thermodynamically stable fluid phase only when the dimensionless density (the 
volume fraction) fj stays below 0.49, with numerical errors reportedly less than 
1%. Numerical simulations^^ also indicate a thermodynamically stable solid 
phase above rj^^ ~ 0.54 all the way up to rj^.^^ = tx^2 l<o ~ 0.7405, the fee crystal 

close packing fraction, although the system may "jam" into a glassy structure.!^ The 

interval (^jV^f^ ) is interpreted as furnishing the coexistence of both fluid and solid 
[941 iV'^l 

phases, i^i^ In the absence of empirical evidence for any further phase transition 
of the hard-sphere system, EH] one may thus assume that the map 7 i— )■ p = ^,(7) 
for a hard-sphere system is a positive, increasing, convex function over M, which is 
asymptotic to a straight line with slope equal to r/^^^ when 7 00. Moreover, the 
map 7 H-)- p = ^,(7) has a kink at 7 = but otherwise is locally real analytic, 
such that away from the kink we have r\ = ^',(7). At 7 = 75^ the left-derivative 
lim^^^j^ pKt) = , and the right-derivative hm^j^^^^ ^',(7) = Vh- Unfortunately, no 
manageable formula is known for the exact ^,(7), but convenient formulas for very 
accurate approximations to ^,(7) are known. 

For the fluid regime 7 G (— oo,7fJ we resort to a formula by Carnahan and 
Starling, E] whose numerological manipulations have led to a simple approximation 
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Pcsil) tc0 P»(7)l-y<7f =• Pmfil) which remarkably accurately fits the empirical data 
obtained in computer simulations. Graphs of the function pes (7) and its derivative 
Pes (7) displayed in Figs.l & 2 of Ref . [33] . 

Definition 2.1: The Carnahan- Starling approximation to ^,^(7) is defined by the 
map 7 I—)- p = pes (7); given by the parameter representation p = g^{fj) and 7 = g2ij]), 



(20) 

where f] is a real parameter in the interval < fj < f]^'^ ^ 0.49. This gives 7fs = 
S-ziVft) ~ 15.208 as right limit for the domain of definition (— 00,7^] o/pes(7)- 
Remark: Note that (fT9l) and ( |20l) are related by a thermodynamic identity for a 
system of many identical hard balls, 

m'M = 9'M- (21) 

Indeed, (|2U]) is obtained from ( IT^ by integrating ( ETj) and conveniently choosing 
the integration constant. As a consequence we have that peg (7) = ^77^ (7) = is a 
dimensionless particle density — as already implied by the stipulated notation. □ 

Remark: Formally (fT9|) and ( 120|) are well defined for all rj G (0, 1), and one may 
want to study this mathematical model in its own right. Whenever we use pes (7) 
as defined for all r/ G (0, 1) we will refer to it as the Carnahan- Starling model, to 
distinguish this mathematical model from the actual hard-sphere physics. □ 

For the solid regime 7 G [7^, 00) we resort to Speedy's effective approximation 
pspil) to p,(7)|^>^j^ =: P,si'l)y whose leading term is determined theoretically while 



^Carnahan and Starling^^ proposed p = g^ir]), with g-^ given in Definition 2.1, as the exphcit 
sum of an approximate virial series for the equation of state p — gt{(jj) for a hard-sphere fluid, 
inspired by the few available terms in the actual virial expansion. Quantitatively their equation of 
state slightly improves over the Percus-YevickSH equation of state under compressibility closure^^ 
(identical to the equation of state obtained from scaled particle theoryS^), from which it differs 
by the extra —rf term in the numerator. 
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the other terms invoke a Fade approximation to fit the numerical simulation dataj^ 
Definition 2.2: The Speedy approximation to p,Xl) ^■^ defined by the map 7 i— p = 

Psp(7); given by the parameter representation p = g.^{T]) and 7 = g^ijf), with^^ 



with a = 0.5921, b = 0.7072, c = 0.601, and 

gM = i.+ r (24) 

where fj ranges in the interval 0.54 ^ fj^ < f] < fj^^^ ^ 0.7402. Note that both g.^{rj) 
and g^ijf) are monotonia increasing on (j]^,rjf^J, diverging t 00 for rj ^ rj'^l^. 

Speedy's paper-^^- features formula ( |23l) . while (IMl) follows from postulating the 



thermodynamic relation rig'^{ri) = g'^{ri) for rj^ < rj < rj^^^; as a consequence, ^',^(7) = 
g'^i'y) for 7 > 7fe. The integration constant is chosen such that 7^ = g^{ri^) ~ 15.208 
is the left limit for the domain of definition [7fs, 00) of pspil)- 
To summarize, we stipulate the following: 

Convention 2.3: In the remainder of the paper, for 7 < 7;^, i.e. in the fluid phase, 
we take ^,((7) := pcs(7) = (S'l ° 5'2^^)(7); '"'^^^ S'l '^'^'^ given by / fTPj) and ^2^) 
in Definition 2.1. For 7 > 7fg, z.e. m the solid phase, we take '■~ Psp(7) = 

(fl's ° fl'7"'^)(7); '"'^^^ S'a ^'^^ gA given by / f^) anc? ( [1^ m Definition 2.2. 

In the next two figures we display the hard-sphere pressure : temperature ra- 
tio (Fig.l) and the hard-sphere chemical potential per particle : temperature ratio 
(Fig.2), both as functions of fj. The second figure in particular will be very helpful 
to consult when reading our proofs in the ensuing sections. 



^Speedy reports that his formula agrees to within less than half the reported error with the 
results of Alder et al.^, which are given by their formula (1), an asymptotic expansion in powers 
of ^/L viz. 

9M = VfL [3 ^„-^^-cp +Ko+ K,{1 - 77/77- ) + 0((^;L - 77)2)] , (22) 
with Kq w —3.44 and Ki w 1 taken from table III in Ref.[3]. 
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Fig.l. Equation of state p vs. 77 for a classical hard-balls continuum. Displayed are fluid branch 
(0 < 7? < 0.49) and solid branch (0.54 < 77 < 0.74) together with the coexistence line (0.49 < 77 < 
0.54; dotted) and the fee crystal close packing limit (broken vertical line at rj = 0.74019). 
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Fig. 2. Chemical-potential-per-particle : temperature ratio 7 vs. volume fraction rj for a classical 
hard-balls continuum. Displayed are fluid branch (0 < 77 < 0.49) and solid branch (0.54 < 7/ < 0.74) 
together with the coexistence line (0.49 < r] < 0.54; dotted) and the fee crystal close packing limit 
(broken vertical line at 77 = 0.74019). 
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We end this subsection by emphasizing that we do need to know what we stip- 
ulated about ^,(7) for the sohd phase even though in this paper we are studying 
only fluid solutions. In particular, "all the hard work" in our paper is caused by 
the following dilemma: to prove a first-order phase transition between two different 
stable fluid solutions one must flnd suitable (0,7) pairs for which has at least 
three all-fluid solutions, but one also must rule out any other globally stable solution 
which takes solid values somewhere in the container. Clearly, a suflicient though not 
necessary condition for the absence of a globally stable solution which takes solid 
values somewhere is the complete absence of any solution taking non-fluid values 
somewhere. This sufficient condition is simpler to implement, but is of course also 
more restrictive. 

The space-uniform van der Waals theory gives a good indication of the difficulties 
ahead. Recall that the space-uniform solutions to the van der Waals problem ([6]) for 
given (a, 7) and ||^||^ < 00 can be graphically determined as the abscissa values of 
the intersection points of the graph displayed in the second figure with the straight 
line ^ I— )■ 7 + a II y II r/. The (0,7) pairs for which a phase transition in the fiuid 
regime occurs while no solid solution exists at all lie in the (0,7) domain which 
corresponds graphically to the family of straight lines rj 7 + a || VH^rJ which have 
three intersections with the fiuid branch but no intersection with the solid branch 
in the second figure Inspection of the second figure makes it obvious that this leaves 
us with only very little "room to wiggle" in (a, 7) space, so that we need to develop 
delicate analytical estimates to accomplish our feat of proving the grand canonical 
gas vs. liquid transition and the petit canonical vapor vs. drop transition within the 
non-uniform van der Waals model for a hard-sphere fiuid with chemical self-potential 
confined to a container. 

3 LOCATING THE FLUID SOLUTIONS IN (a, 7) HALF 

SPACE 

In this section we give some sufficient and some necessary a priori conditions con- 
cerning the existence of solutions 77 of ([8]) which do not take values outside the 
fiuid regime, i.e. for which 7 — {aV * 77)^ < 7^. We shall write V * t] for either 
{V * 7])^ or {y * if)^^ whenever it is clear from the context what we mean. We set 
\\V * l||c{;'(A) = ^a; where || ■ ||cO(A) denotes the uniform (supremum) norm for C°(A); 
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notice that ^3 = We also introduce the notation = {t] : \\ri\\ < ^} for 

the closed ball of radius ^ in C°(A). 

We begin with some sufficient conditions for existence. 

Proposition 3.1: Assume that the inequality 

-f + a^^r]-gM<0 (25) 
is satisfied for at least one fj G (0,r/j^], so that the algebraic fixed point equation 

r] = g;\l + a^AV) (26) 

has at least one solution G (0, fj^'^] . Let r/™ be the minimal and Tjf the maximal solu- 
tion in [0, r/j^] of / f^) . Then in the truncated positive cone C°^^(A) HiS^M there exists 
a pointwise minimal and a pointwise maximal fluid solution of denoted r]'^{'r) 
andr]^M{r), respectively. In particular, the iteration sequences {i]^"^}'^^=o defined by 

^(n+i) _ g-i^^ -aV* r/(")) (27) 

with starting densities r]^^^ = r/™ and t]^^'^ = g~^{'j), respectively, both converge point- 
wise to the minimal solution rj'^ir), the former monotone downward and the latter 
monotone upward. Starting the iteration map ( [F7| j with r]^^'^ = r/*^ yields a sequence 
which converges pointwise monotone downward to the maximal solution r7^1,j(r). 

Remarks: (a) Since, by hypothesis, fl25l) is satisfied, and since g^ is continuous with 
lim^j^o (^) = —00, the straight line 1— )■ 7 + a^^f] intersects the curve f] 1— )■ g2(j]) 
at least once (and at most three times) in (0, 77^^]. Therefore a maximal point of 
intersection r/" < r/^^ does exist, (b) Proposition 3.1 does not state that r^^M is the 

maximal solution in C^^ fl i3^<; however, 77^ is automatically the minimal solution 
in C^_|_ n Bj^<. (c) Maximal and minimal solution, tj^'m and 17™, may coincide. □ 

To prove Proposition 3.1, all we need to know about \^ is ^ G -^^^(A) and < 0. 

Proof of Proposition 3.1: Consider first the case 77^'^) = fj^^. Since 7 g~^{'y) is 
strictly monotonic increasing, and since —{V * l)(r) < \\V * 1|| = but with 

-{y * l)(r) ^ <1)^, the iteration (E?]) yields r/(")(r) < ri^'^-^\r)\/ n G N, and even 
r/("+i)(r) < r7(")(r) Vn G N and all r G A. Since 

g;\l-aV*r^)>g;\-l)>Q, (28) 
15 



the iterates are bounded below by a positive number. Hence, the iterates converge 
pointwise down to a strictly positive function tiI'm which clearly is entirely fluid. By 

'''a 

the C°(A) continuity of the operator g^^{;~i — aV * ■), the function rjlii solves ([8]). 

As in Ref.[6] it can be shown that ?7*m is the pointwise maximal solution in 
C°^^ n Bj^M. For suppose that r] < fj" is any solution of (|H]), then by the monotonic 
increase of g^'^il ~ ciV * ■ ) and by the fact that r/^^ is a strict supersolution for (|8]), 
we can conclude that gT^^i^l — cxV * ■ ) maps \r],rf-l\ into itself. Therefore, r] < rj^ii, 

and this proves that ?7*m is the pointwise maximal solution in C°_,_ fl B^jm . 

By essentially the same arguments, starting the iteration with 77*^°^ = rj™ yields 
a monotone downward converging sequence of iterates with limit r/™, and 17™ is the 
pointwise maximal solution in C^_,_ fl Bjf^. 

Next consider the case r]^^^ = g~^{'y)- Using again the strict monotonic increase 
of 7 I—)- (7~^(7), this time combined with the positivity of —{V * l)(r), we conclude 
that the sequence (1271) iterates pointwise monotone upward. By fl2Sl) all iterates are 
strictly positive. Moreover, by induction it follows that, if i]'-"^ < r/", then 

= g;\7 - c^V * < g;\^ + = C • (29) 

Clearly, r]^^'> < rj'^^ so the sequence is bounded above by fj^. It now follows that 
it converges pointwise to a strictly positive solution ^iX" — ''^a of dHD, and also that 
this solution is entirely fluid. Moreover, similarly as for the maximal solution it now 
follows that ?7™ is the pointwise minimal solution in C°_,_ fl -B^j™, hence in C^_^ fl i3^<. 

Lastly, the proof that ?7™ = r/" is a minor variation on the proof of Corollary 4.5 
in section IV. ■ 

By a slight sharpening of ( l25ll we can improve Proposition 3.1 to the following. 
Proposition 3.2: Assume that 

7-7, + a<l>^r/< <0. (30) 

Then ( fl3]) is satisfied for fj = rj^, so Proposition 3.1 applies. Now the pointwise 
maximal fluid solution rj^ii of ^ in C^^ fl B^^m is in fact the pointwise maximal 

fluid solution m C°_^ fl i3^< . 

Proof of Proposition 3. 2: Since 7^ = f/al^fe )' (00]) implies that f l25|) is satisfied by 
fj = f]^, so all conclusions in Proposition 3.1 apply. 
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To show that the pointwise maximal solution in C°_,_ni3^M is in fact the pointwise 
maximal solution in H Bj^< , we notice that rj < is a strict supersolution a.e. for 

([8]). This implies that the sequence {r]^"'^^^Q defined by (1271) with initial value 
7^(0) = f^< iterates pointwise monotonically downward, strictly monotonically a.e., to 
the pointwise maximal solution in C°^^ fl Bjj< of ([8]). We show that this solution 
is in C^_^ n B^M, and so, a forteriori, it is also the pointwise maximal solution in 

By fl3UI) . T^j^ is a supersolution for ( l26i) . Therefore, either r/^^ is itself the largest 
fixed point in [0,r/j^] of fl26|) . or else the sequence 

}^=o defined by 
¥"+'^=C'(7 + «<^'a¥"^) (31) 

with initial value r]^^'' = rj^^ iterates strictly monotonically downward to the largest 
fixed point in [0,rjj^) of fl26l) . which in either case is f]-<. Moreover, with f]^'^^ = r/^^ = 

^fs 

7^(0)^ for each n > we have 

r/(") < rj(") , (32) 
because r/*^"") < rj^"") for some uq >0 implies that 

^(no+i) _ -aV* r/("o)) < g;\-f + a^^r^^"")) = r/("«+^) . (33) 
We conclude that 

ri^'< := lim r/^") < lim rj^") = rj;^' . (34) 

Hence, r]^< = tj^m, so ?7*m is the pointwise maximal solution in fl Bj^<. ■ 

Remark: For our V, the dominance can be sharpened from "<" to "< a.e." by 
noting that obviously r]^^^ < r]^^^ a.e. □ 

If we consider the extension of ([8]) to all 7 G M, with p',( • ) = pcg( • ) = 5'^^( ■ ) 
for ■ < with p'^ meaning left derivative, and with p',( ■ ) = g^^{ ■ ) when ■ > 7^3, 
with p', now meaning right derivative, covering fluid and solid branch as explained 
in Convention 2.3, then the existence results of Propositions 3.1 and 3.2 can be 
complemented by a result about the non-existence of solutions to the so extended 
([8]) which are not fluid somewhere in A. 

Proposition 3.3: // the inequality 

7-7. + «$^r7;L<0, (35) 
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holds, then the extended fixed point problem Ij^ does not have any solution that takes 
values outside the hard-sphere fluid regime somewhere in A. 

Proof of Proposition 3.3: Since rj < rj'jl^, and since g^iji) > lis for all rj G (j]hyVfcc]y 
we conclude that ( 135!) implies for all r/ G (fj^ , rJ^^J that 

7 + a<l>^f7 < ^4(r?). (36) 

Now suppose a solution rj of the extended ([8]) would exist which in some open 
subdomain of A is solid. Then, clearly, r}^ < ||?7||c'0(a) ^ ^/L' ^^^1 since 7 1— g^'^i'y) 
is increasing, we conclude that in the solid region (i.e., in A^) we have 

ll^ll,o<,,< ^4-^(7 + «$^ 11^11,0, J (37) 



as a consequence of the extended ([8]). But by applying g^ to both sides of fl37|) . this 
leads to a contradition with f p6|) . Hence, no solution of the extended ([8]) can exist 
which somewhere in A is not a hard-sphere fluid. ■ 

The next result requires V G L^(R^). It relates the algebraic fixed point problem 
for constant solutions in M.^ of ([7]) to the problem ([8]) in bounded A C M^. 

Proposition 3.4: Let a \\V\\^ = A^Itt"^ /Ax^) + AY(4'7r/K^). Suppose the algebraic 
fixed point problem ^ has a solution r/^^^ < r/<, so that fj^^^ satisfies 

r) = g-\^ + a\\V\\^r)). (38) 

Then for all domains A C R'^ the fixed point problem ^ with aV = AwVw + AyVy 
and p,{'~f) given in Definition 2.1 has a hard-sphere fluid solution. 

Proof of Proposition 3.4: By subadditivity of the norm, we have 

<Aw||Vw*l|| +Ay||V'y*1|| • (39) 

C0(A) C0(A) eg (A) 

Since V^{\ ■ |) G ^^(R^) and Vy{\ ■ |) G L\R^), we have 

||K.*1|| „ < ||Kv*l|L 3 = (40) 

in 

\\Vy*l\\ < \\Vy*l\\ (41) 

C0{A) c0(k3) ,^2 
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With dSni), din]), and dH]), we thus have 



a ||r * l|l < A^-— + A^— = \\Vl , (42) 

vahd for all A C M^. Hence, if (138|) has a solution rj^^^ < rjf'^, then by ( l42ll this r/^^^ 
is a super solution for and Proposition 3.1 now concludes the proof. ■ 

We turn to the necessary conditions for the existence of fluid solutions. 

Proposition 3.5: // the inequality 

7-7, + a$^p:(7) >0, (43) 

holds, then the extended ^ does not have a hard-sphere fluid solution. 

Proof of Proposition 3.5: Since ^ < and a > 0, and since ^^(7) > 0, it follows 
directly from ([8]) that any solution 77 of the extended ([8]) satisfies the lower estimate 

r^iv) > p:(7) (44) 

for all r G A. Convoluting (jHj) with —V (> 0) and multiplying by a yields 

- {aV *ri){r) > -{aV *l){r)p',{-f) (45) 
for all r G A, from which it follows that 

7 + \\aV * vWc^iA) > 7 + a^Ap'.il)- (46) 

If (H3|) holds, then from P6l) it follows that 7 + ||«V^*?7||cO(A) > 7fs, and so ||?7||cO(A) > 
f]^. Therefore, violation of (143|) is a necessary condition for the existence of an all 
fluid solution of ([8]). ■ 

We conclude this section with an obvious non-existence result. 
Proposition 2.6: // the inequality 

7-7f. >0, (47) 

holds, then the extended ^ does not have a solution which is fluid somewhere in A. 

Proof of Proposition 2.6: Trivial. ■ 
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4 A (a, 7) REGION WITH UNIQUE FLUID SOLUTIONS 

We now locate a connected region in (a, 7) space in which there exists a unique fluid 
solution for each pair of {a, 7) parameter values. The pertinent unique fluid solution 
need not be the unique solution per se, yet any other solution of ([H]) would necessarily 
take nonfiuid values somewhere in A. 

Our existence and uniqueness results are based on the following theorem, for 
which much less is assumed about ^,(7) than stipulated in Convention 2.3. 

Lemma 4.1 Consider ^ for a map 7 1— )■ ^,(7) of class 00, 7) which is strictly 
positive, increasing, and convex, and for which 

K{^):= sup p:(7)<oo. (48) 

7e(-oo,7) 

Assume 7 (< 7) and a{> 0) are such that the operator p'^i^j — aV * ■ ) maps C^^nBfj 
into itself, where Bfj = {rj : \\ri\\ < fj] is a hall of radius fj = ^',(7). Assume 

furthermore that 

i^(7)a<f ^ < 1 , (49) 
with = \\V * 1\\ ^ , as defined above Proposition 3.1. Then there exists a unique 



solution G C°_,_ n of In particular, the iteration sequence starting with 



any 77*^°^ G C^_^_ D Bfj, converges strongly in C°(A) to the unique solution. 

Remark: Lemma 4.1 improves over Theorem 6.4 of Ref. [33], where uniqueness and 
strong L} convergence are established under the same condition ( l49l) . □ 

Proof of Lemma 4-1- By hypothesis, the operator p',(7 — aV * ■) maps the 
II ■ llco(A) closed set C°^_ fl Bf, into itself. This implies that 7 — aV * 77 < 7 for any 
7/ G C^j^ n Bf^. This together with (HH!) in turn implies that p^'(7 — OiV * rf) < K{'y) 
for any 7] G C^_+ (iBfj. 

Consider now two sequences {r/J"^ G nBr^}^^o,t = l,2, defined by ((271), with 

rj^^ 7^ ^72°^ on a fat set. Set —{V * ?7-"^)(r) = ^^■"■'(r). Pick any 1 < g < cxd. Then, by 
the fact that ^','(7 — aV * rj) < K{'y) for any rj G C^_^_ fl Bfj, we estimate 



(n+l) (n+1) 

m -vi 



<? 


L 







p:(7 + «#(r)) -pI (7 + «0j"Hr) 



1 „ 



20 




4"'{r) 




dip 



(Tr 



4"' 



L9{A) 



(50) 



By the definition of (t)f'\ followed by an obvious estimate and then by an application 



of Holder's inequality with conjugate exponents g, g', we estimate 

j |0^")(r) -0S")(r)|''dV = j j -V{\T-i\) (r]f'^-r]t^'^ (i) d'f 



< 



< 



-vy' * 1 



g/q' 



Lt/l' (A) 



(n) (n) 

V2 -m 



(n) (n) 
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LI (A) 



(51) 



Combining ( 15U]1 and gives, after taking the gth root, 



(n+l) (n+1) 




i-vy' * 1 




(n) (n) 






L1{A) 




Ll/l'(A) 




LI (A) 



(52) 



for all q G (1, oo). By taking g — )• oo, and noting that here ess sup = sup, we get 



(n+l) (n+l) 



< Kh) \\aV*l\\ 

n ~ C0{A) 



(n) (n) 

Vl -V2 



(53) 



eg (A) 



By hypothesis (H^ . we have K{'y)a ||y * 1|| < 1, whence from (153 p we con- 

C|,'(A) 

elude that the map t] p',(7 — ctV * r/) is a C° contraction map in the closed 
truncated cone fl Bfj. We now apply the contraction mapping principle^^M 
and conclude that a unique fixed point of rj ^',(7 — aV * rf) exists in C°_,_ n Ef^. In 
addition, the proof of the contraction mapping principle implies the C° convergence 
of the iteration sequence ([27]) for any initial density rj'^'^^ G C°_^ H Bf,- ■ 

We now return to our ^,(7) given by Convention 2.3. In our first application of 
Lemma 4.1 we set 7 = 7fs (~ 15.208). The following input from Ref.[33j capitalizes on 
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the fact that the graph of 77 (^) has a unique inflection point at 77 = ?7; ~ 0.130. 
Lemma 4.2: The regular global maximum K{'~fis) of ^"3(7) over the set (— 00,75^) 

occurs at 7; ^ —0.67 ai which rj^ = gf^"^ (7;) 0.130 and p'csili) — ^{its) ~ 0.047. 

We are now in the position to state the following Corollary of Lemma 4.1. 

Corollary 4.3: Let the parameters (a,7) satisfy the hound l[25\) . and let a satisfy 
the inequality \\aV * -'-II ^ 21.20. Then there exists a unique fluid solution of 

Proof: By hypothesis, the parameters (0,7) satisfy fl2S]) . This implies that the 
operator ^^3(7 — aV * ■ ) maps fl i3^< into itself. Next, using Lemma 4.2 and 
1/0.047 ^ 21.20, we conclude that' ||a1/*'l|| ^ < 21.20 implies g9]). Lemma 4.1 

now guarantees us a unique solution G fl i3^< of (JH]). ■ 

It is interesting to compare (1491) to the sharp criterion for uniqueness, irrespective 
of 7, of a solution rj < rj^ to the associated algebraic fixed point problem (126|) . 
Geometrically, this criterion for uniqueness is that the slope of the straight line 
77 I— !■ 7 + '^o^ surpass the smallest derivative of the curve fj h-> g^ijf)^ or 

< g'^iVi)^ with rjj ^ 0.130 defined in Lemma 4.2. From the definition of pcsil) 
we then see that this criterion is precisely i^(7fja$^ < 1, with ii'(7fj = Pcsili) (~ 
0.047) given in Lemma 4.2. Thus, f l49|) is the direct analog of the geometric criterion 
that governs the associated algebraic fixed point problem fl26|) . except for the case 
of equality i^'(7fs)a$^ = 1, about which the contraction mapping principle is silent. 

If K{'jf^)a^^ > 1, then there exist values of 7 for which fl2Bl) has either two or 
three solutions. In that case we can still arrive at a uniqueness theorem for fl2Bl) under 
the condition on 7 that it be not too large. Similarly, if (H9|) is violated. Lemma 4.1 
still gives a uniqueness result for ([8]) by appropriately restricting 7 from above. For 
this second application of our Lemma 4.1 we introduce the following. 

Definition 4.4: Given A, for each a we define 7a (a) to be the largest upper bound 
on 7 such that for each 7 < 7^(0;) there exists a unique positive solution 77(0;, 7) of 

Remarks: (a) Since g',^{vi) > and S'allO^^fe]) = ("00,7^], clearly 7^ > —00; (b) 
7A(a) has a discontinuity when a$^ii"(7(J = 1. □ 

Corollary 4.5: Let a satisfy i^'(7fja;$^ > 1, and let 7 < 7A(a;). Then fj{a,'y) < rj^, 
and ^ has a unique fluid solution rjjy G fl Bjj< ; in fact, rj^ G fl B-i^aa) ■ 
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Moreover, the iteration sequence defined by (7 — aV * 1]^"^), starting with 

any i]^^^ E C°_^_ fl Bjj<, converges in supnorm to this unique fixed point. 

Proof: Since ^(7^)0;$^ > 1, by definition of 7a we see tliat fj{a, 7) < r/,. Tlierefore 
all fj G supersolutions for fl26|) . and thus strict supersolutions for ([8]). By 

the type of argument presented in the proof of Proposition 3.2 we conclude that no 
fluid solution of ([H]) exists which is somewhere larger than r/(a,7). 

Now pick any rj G C^^ fl Brj(a,-y)- Since 7 < 7a(«), (fi'7^)'(7) > 0, < 0, we have 

< ^7^(7 + a$^fj(a, 7)) =77(0,7). (54) 

Therefore, the operator g^^ij ~ * ' ) maps fl Bjj(^a,'y) into itself. 

We observe that g^'(j]{a, 7)) > so that f/(a, 7) is a stable fixed point of ( !26|) . 
The stability of 7/(0,7) ^ind the convexity of g~^{i^) for z/ < 7 + a$^r/(a,7) implies 
that 

ir(7A)a$A < 1 , (55) 

where 

K{%) := supj {g~^y{-f + z^) : 7 G ( - 00, 7A(a)) A u < a$^rj(a, 7) } (56) 
We now can apply Lemma 4.1 to 77 G C^_^_ fl Brj[a,'y)- The proof is complete. ■ 



5 A {a, 7) REGION WITH SEVERAL FLUID 

SOLUTIONS 

When V G L^(M^), it is readily shown that there is a connected region in (0,7) 
parameter space in which the van der Waals' algebraic fixed point equation ([6]) for 
constant density functions in has three solutions inside the fluid regime, fj"^^^ < 
^vdw < ^vdw — Vity so these solutions satisfy ( l38l) . The smallest and the largest 
ones are stable under iterations while the intermediate one is unstable. Intuitively 
one expects that when A C is a container of macroscopic proportions, and k,~^ 
and >c~^ are molecular distances, then for most (0,7) in the three fluid solutions 
region for the algebraic f l55]) our nonlinear integral equation (|H]) should also have a 
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small and a large fluid solution which are stable under iterations, while the unstable 
solution rj"^^ of fl38|) should be replaced by an interface type solution of ([8]) which 
is unstable under iterations. Numerical integrations of ([8]) with \^ = Vw for a ball 
domain A = S/j with moderately large R = 50/x support this expectation. A 
rigorous proof is desirable. 

In this section we use monotone iteration techniques with sub- and supersolu- 
tions to show that at least part of this multiplicity region for the algebraic equation 
(l38l) corresponds to a multiplicity region of the integral equation ([8]) — for certain 
sufficiently small A. We will prove that at least three hard-sphere fluid solutions 
exist in some region of (a, 7) parameter space, two of them stable under iteration 
and one unstable. We will not show that exactly three fluid solutions exist; in fact, 
it might not be true that exactly three fluid solutions of ([S]) exist whenever it has at 
least three fluid solutions. 

Recall that the starting function r/*^*^^ = g~^{'y) is a subsolution for in any A, 
and it launches an iteration sequence which converges upward toward the pointwise 
minimal solution; see Proposition 3.1. We also know from Proposition 3.4 that 
when r/^j^ < ry^^, then any starting function rj^'^^ G [ntdw^Vk] is a supersolution 
for ([H]) in any A, and it launches an iteration sequence which converges downward 
toward the pointwise maximal fluid solution. One can rule out that the pointwise 
maximal solution coincides with the pointwise minimal solution if a sufficiently large 
subsolution of ([8]) in A is available from which the iteration ry("+^) = p'^^j — aV *r]^''^^) 
converges upward toward a fluid solution which is larger than the pointwise minimal 
solution to ([H]). 

Constructing suitable subsolutions that imply a [a, 7) region of multiple hard- 
sphere fluid solutions is a very difficult business, yet much easier for the Carnahan- 
Starling model. We will take advantage of this fact and, until further notice, first 
discuss ([8]) with §3,(7) replaced by pes (7) for all 7 G M, viz. 

v{r) = p'as{l-ic^V*r^Ur)). (57) 

Subsequently we seek those solutions which nowhere in A surpass r]^ . We emphasize 
that our multiplicity results for the Carnahan-Starling model in general have no 
bearing on the hard-sphere fluid; however, there will be a small sliver in (a, 7) space 
for which our Carnahan-Starling multiplicity results yield multiple hard-sphere fluid 
solutions. 
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So recall that ( ■ ) = ( ' ) consider the algebraic fixed point problem 

V = 97^il + OiTV) (58) 

for 7 G M and a G M+, where r G M+. Multiplicity of solutions of ( 158|) can only 
occur if a is large enough, namely (recalling Lemma 4.2 and Corollary 4.3) if 

ar > min g'Cf]) = g'Cf],) ^ 21.20. (59) 
77e(o,i) 

In addition, 7 needs to satisfy 7 G (7cs (o^'^)' 7cs ('^''"))! with upper and lower interval 
Hmits given by 

7cl("^) = 92{r]<) - arr]^, (60) 
7cl("^) = 92{V>) - a^^>, (61) 
where < are the two distinct solutions to the equation 

ar = ^;(rj), (62) 

which exist only when f l59|) is satisfied, in which case r/^ < r/j is a decreasing, and 
r/^ > rjj an increasing function of ar. While it does not seem feasible to write down 
closed form expressions of the functions ar h- )■ f] ^ and ar 1— )• 77^, their asymptotics 
for ar ~ g'^ijii) (recall ( l59l) ) and ar ^ Q^^i) can easily be worked out, which gives 
us 



7cl(«^) 



92{Vi) -Vi^r 
-ln(ar) - 1 + 0[l/ar] 



7S|(ar)x -^^"^ -^^"^ .n^.) 

-far + 0([ar]3/4) 



ar > ^^(rj,) 
ar > g'^{%] 



(63) 
(64) 



where we used the identity g[{j]^ = r]^g'^{ri^) to simplify. Numerically, g'^'{rj^ ^ 
1235.22. 

So the algebraic fixed point equation ( l58ll has three solutions for all (a, 7) G 
eg(r), where eg(r) = {(a, 7) : ar > g',{rii) ^ Tcliar) < 7 < 7cl(ar)}. Note 
that the boundary 5B^'|(r) is given by two functions of a which depend on a only 
through the product ar. Hence, for fl58|) . triple solution regions in the (a, 7) half 
plane for any two different r = ri and r = T2 differ from each other only by some 
scaling along the a axis, viz. they are affine similar. Since for fixed r the upper 
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boundary curve 7^'|(ar) diverges to — oo logarithmically while the lower boundary 
curve 7^'|(ar) does so linearly when a becomes large, it follows that any pair of triple 
domains Ocs(''"i) ^^^d Q%1{t2) has a non-empty intersection. 

We next identify functionals of V which can be substituted for r to construct 
sub- and supersolutions for f lFTI) . Since both our van der Waals kernel Vwdi"!) and 
the Yukawa kernel Vydrl) are monotonic increasing negative functions of |r| = r, for 
aV = A^Vw + AyVy and any container A with diameter 0(A), we have 

V{\r -r'\) <V{0{A)) V r, r G A. (65) 

We define the abbreviation 

VI/, :=-V(0(A))|A|. (66) 

Subsolutions for (1571) can be constructed by setting r = supersolutions by setting 
r = or r = Note that for bounded A C R'^ we have the chain of inequalities 

VI/, <$^3 = ||l-||^. (67) 

Since our findings about the triple algebraic solutions domain for (l58l) imply in 
particular that for any bounded domain A C we have 6^'|($,) fl 9|^'|(vE',) 7^ 
and also ©csdl^lli) ^ ©cs(^a) '^^^ show, with the help of monotone 

iterations and the mountain pass lemma, that for each (a, 7) G 6^'|($,) fl 0^'|(v&,) 
and each (a, 7) G 0cs(II^IIi)'~'®cs(^a) ^^e fixed point equation (|57|) has at least three 
solutions in C°(A), which are ordered. However, for physically interesting domains 
A the sets OcI('^'a) ^ ©cII^a) and Ocldl^lli) H 6^'|(^,) are generally very bad 
approximations to the full set of such (a, 7) points. The reason is that for physically 
interesting, i.e. macroscopic domains A, the ratio ^1/,/$, is tiny, converging to zero 
as A t R^- Worse, 6p|($,) fl 0^'|(vEf,) may even be a totally useless estimate of the 
three hard-sphere fluid solutions regime of ([8]), in the sense that the largest solution 
of (1571) obtained by this method may always take values outside the physical range 
of hard-sphere fluid densities. 

The following variation on our strategy yields more desirable multiplicity results. 
For bounded A C M^, let <;A C A denote a rescaling of A into A by a factor ^ < 1, 
so that 0(^A) = <;0(A) and |^A| = ^3|A|. Then for aV = A^V^ + A^Vy the map 
<^ H> vp^^ = — ^^y((; 0(A))|A| takes a global maximum at ? = <^ (which might not 
be unique; it is unique when V = V^^f or V = Vy)- We always mean the largest 
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Suppose now that A is a container domain of macroscopic proportions, and that k ^ 
and are molecular distances. Then <j ^ 1, and \/ q > ^ we have the ordering 

For spherical macroscopic A (see Appendix A), and with n = h = 1/2, we have 
~ 50\E'^^, while the first inequality separates two quantities "a universe apart." 

Proposition 5.1: Let aV = A^Vyf + AyVy- Let A G M.^ be a container for which 
tK C A. Then for each ^ G [<j, 1] and (0,7) E ©cII'^'.a) ^ ©cII^^a) equation 

v{r) = p'csil-i(^V*vXd^)) (69) 

has at least three distinct solutions in C^(^A). In particular, has a pointwise 
minimal and a pointwise maximal solution, both of which are stable under iteration, 
and a third, unstable solution which is sandwiched inbetween. 

Proof: For each (0,7) G ©csI'^s-a) ^ ©cs(^4a) ^^e algebraic fixed point equation 
( l58l) has three solutions for r = and for r = \1'^^, denoted r/^ < ff^^ < rj" 
and r/™ < rf^^ < r]fj^, respectively (suppressing their dependence on (a, 7) from being 
displayed). Moreover, since "^.^ < we have rj!^ < rj^ and Tjf^ < TJ^^; the ordering 
of the unstable solutions is rf^^ > rf^^, but this is irrelevant for our arguments. 

Now consider the iteration r/("-+^) = p'^g^'j — (aV * ^^"'^Xa) ^fe,+('^^) ^ with 

i°^(r) = Pcs(7-c£«^(|r-f|)^'^~) VrG^A (70) 

and either /x = m or M. It is easily verified that rfif'iv) is a subsolution of ( l69l) . Since 
Pcsi ■ ) > p'(.g( ■ ) < 1, each rif\Y) launches a monotonic increasing sequence 
{?7^"'''}^g G C°^^(A) n i3i which converges pointwise to some solution ?7^(r) of ( l69l) . 

Moreover, we have (r) < fj^ir). To verify this claim, we note on the one hand 
that in Proposition 3.4 we already showed that the constant function r i— )• rj™ is a 
supersolution of (!69|) for any ^A (recall, this follows from —{V* 1)^ < $^ for any A), 
so that with r/™ < fj'^^ we find r]m{r) < r/™ , and now we conclude as in the proof of 
Proposition 3.1 that ?7^(r) < ^J^; incidentally, < rjj. On the other hand, '?7M''(r) > 
^«A > ^7;Vr G ^A, and since the iteration {r^M^j^o monotone upwards, it follows 
that vlti''^) > ''^K^)^'^ ^ addition, r/i/^r) > r7m''(r)Vr G ^A, so the strict 
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monotonic increase of the iterations now guarantees that '7^^(i') > ''7l(r)Vr G <;A, 
and since '7^^(r) > //^(r) Vr G ^A, it even follows that ?7^^(r) > //^(r) Vr G ^"A. 

Standard results about monotone iterations in ordered Banach spaces show that 
r/^(r) and '7^^(r) are stable under iterations, and also locally V stable; see Propo- 
sition 3.1 in Ref.[B]. The existence of a third, unstable (under iterations and in V 
sense) solution sandwiched between 77^(r) and vlii^) follows, via the mountain 
pass lemma, from the local V stability of ^7^(r) and rj'^ (r) and the strong C^(A) 
differentiability of the functional V^^^lr]]. 

Lastly, a forteriori the unstable solution sandwiched between 77^ (r) and rj^^ (r) is 
also sandwiched between the pointwise smallest and the pointwise largest solutions, 
?7™ (r) and 'q^Hr), of ( l69l) . obtained by the iteration r]^^+^^ = ^^3(7 — (aV * V^"^Xa) 
from, respectively, 77^°^ = p'csil) any r]^^^ = fj^^^ > f]^[] cf.. Proposition 3.1 with 
(0,r7<] replaced by (0, r/^^] or by (0, 1)), and which are stable under iterations. EI ^ 

■ 

Our proof of Proposition 5.1 reveals the ordering 

Pcs(7) < < <(r) < ^^^.(r) < <(r) < tl- (71) 

Our next proposition shows that the first "<" actually is an identity. 
Proposition 5.2: Under the hypotheses of Proposition 5.1, we have 

V:Ar) = Vlir) (72) 
Proof of Proposition 5.2: An obvious variation on the proof of Corollary 4.5. 

■ 

Remark: We surmise that also ?7^(r) = rj^'/ir) but have not been able to prove it. 

□ 

For a macroscopic container A, Proposition 5.1 tells us in particular that the 
Carnahan-Starling model (169|) with <j = 1 has at least three ordered solutions when 
(a, 7) G Ocs(^a) ^ ^cli^ih)- C)ne of these solutions is bounded above by while 
another one takes (some) values larger than r/,. For large enough a and negatively 
large enough 7 (recall that 6^g($^) n6^g(\E'^^) is unbounded) this large solution will 
take values larger than r/^^, possibly even larger than r/^^^ ^ 0.7402. Those solutions 
do not seem to have an interpretation in terms of hard-sphere systems. 

We now return to our task of finding multiple solutions of ([8]) which all take 
only hard-sphere fluid density values. Unfortunately our analytical control is not 
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good enough to find a subset of 6j^|(<l'^) fl 6^'|(\E'^^) which satisfies our wishes, and 
it's even more hopeless to naively seek an admissible subset of 9^'|($^) fl 6^s(^E'a)- 
However, if we shrink the size of A by choosing a suitable <^ G (^, 1), then we can find 
a subset of 6cs('^<;a) I^0cs(^5a) which at least three solutions of fl69|) take values 
only in the hard-sphere fiuid regime, i.e. for which rjf < r/^^ ^ 0.49. So we impose 
the restriction fj^^ < rj^ on 6cs($^a) ^ QcsI^^a) ^^"^ ^^^^ admissible 

To analyze the effect of this restriction we impose it on QqI(t). Let (r) denote 
the (a, 7) domain featuring three solutions of (l58l) in the hard-sphere fluid regime. 
Recalling the proof of Proposition 3.2, it is readily verified that Q^fi'^) is given by 
e^'f (r) = {(«,7) : g'^{r),) < ar < g'^{r]<) A 7.t(ar) < 7 < 7.t(ar)}, where 

7.t(«^)= Tclic^r), (73) 

7.t(«r) = min{7S'|(ar) , 7^^(«^)} , (74) 

with 

7ff(ar)=7f.-^^«r. (75) 

We note that the two boundary curves jtHar) and jll^lar) intersect at the endpoints 
of the allowed ar interval, i.e. at ar = g'^(j]^) and ar = g'^iv^)- ^Iso the boundary 
90^'/^ (r) is given by two functions of a which depend on a only through the product 
ar, and this implies for (!58|) that also triple hard-sphere fiuid solution regions in the 
(a, 7) half plane for any two different r = ri and t = T2 differ from each other only 
by some scaling along the a axis, i.e. once again these triple regions are affine similar. 
However, distinct from the set 6^'|(r), the set (r) is bounded, and since it is also 
bounded away from ar = 0, if ri and T2 differ by too much then (ri)n6^'f (r2) = 0. 

Thus, to carry out our construction of subsolutions presented in the proof of 
Proposition 5.1 we need to limit the size of ^"A to make sure that ©^'f (^^a)/®»?('^5a) 
is not too small. Since A is supposed to be a macroscopic container domain, this 
means that <j > <^ has to be chosen sufficiently small. Recall that the maximum of ^'^^ 
then occurs for one or more ^ <^ 1, and we stipulated that we mean the largest ^ in 
case ^ is not unique. We can precisely, though only implicitly characterize the range of 
scaled domains <jA for which our construction of subsolutions presented in the proof of 
Proposition 5.1 can be carried out. Namely, the intersection ©^'f (^^a) '"'^•/^('^^a) ^ 
for all q G [^, <^), where <^ > ^ is the unique <^ value for which the lower boundary of 
©•/^(^^a) o^iy touches the upper boundary of ©^'/^($^a) (possibly more than once), 
determined by 

7.t(avl>^J = 7.t(a<|.^J , (76) 
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d^^:na^J= 9,7.?(«$^J. (77) 

The upshot is: 

Proposition 5.3: Let ^ G [<?,<") and (0,7) G Ql'li'^^J n 6^'f(<l>.^). T/ien 

^(r) = p:(7-(«^*^),A(r)) (78) 

has at least three ordered solutions in C°(^A) fl -B^<, two of which can be computed 
by iterating with r.h.s.^^7^, starting from ([7^ with fi = m and fi = M , respectively. 
Remark: It is helpful to have a geometric illustration of the situation. Recall that 
9^j'^(r) is the bounded domain in (a, 7) half space determined by (!73|) , fl74l) . (175!) for 
which the algebraic fixed point equation (158!) has exactly three solutions in the hard- 
sphere fluid regime. For the various r > values associated with A which we have 
encountered in this section, all the domains Q''^1{t) are located in the negative 7 half 
of (a, 7) half space. They have roughly the shape of a receding moon crescent, being 
affine similar to each other by horizontal scaling (along the a axis). The domains we 
have encountered are arranged as follows: 0^'/(||^|| J is the leftmost domain, followed 
by e^j'^($^), then ($^.a), then e^f'^($.^), and finally e^j'^(^,^). For macroscopic 
A we have Qtli^f.) ^ (ll^lli) 7^ 0. in fact 0^'/($a) ~ (ll^lli). and we have 
(^.a) ^ 0^?(*iA) 7^ for all q e [^,4)] however, H QlH^^J = 0, and 

there is much space inbetween. □ 

For general macroscopic domains A it is not easy to come up with good explicit 
estimates on but in our section on spherical domains we will see that <fA is not 
exactly what one would call a macroscopic domain. So Proposition 5.3 falls far 
short of our ideal goal, which is to construct suitable subsolutions in macroscopic 
A which imply that for most if not all (0,7) G 0^j*^($^) equation ([8]) has (at least) 
three solutions whose range is in (0,r/j^). On the other hand, with the help of 
variational arguments we will be able to show that for a significant fraction of pairs 
(a, 7) G ($a) fixed point equation ([H]) has at least three solutions whose range 
is in (0,fjf^), indeed. These arguments invoke our functional Va^^[r]] given in ([9]). 

6 V STABILITY AND THE GAS o LIQUID PHASE 

TRANSITION 

Consider first V G L^(M^) and recall that 0^'j*^(||V^|| j) is the bounded domain in [a, 7) 
half space determined by ( 173|) . ( 17il) . ( I75|) with t = ||V^|| for which the algebraic 



30 



fixed point equation (jS]) has exactly three solutions in the hard-sphere fluid regime 
which are spatially uniform solutions of ([7]). This triplicity region of uniform hard- 
sphere fluid solutions contains a phase transition curve 7 = Jlf^ia) along which the 
mean pressure functional Ila,-y{ri) := limA^-RS |A|~^'P^,^[r7] has an uncountable family 
of global maximizers for each (a, 7) = {a,'j'^f^{a)) — the variational problem for 
^a,-y{v) is degenerate! Amongst its global maximizers are a small (^J;™^) and a large 
(^vdw) spatially uniform solution of ([7]). For spatially uniform density functions fj, 
the functional Ha^'y takes the simple form 

= P.(7 + a lll^ll, rj) - i« II f , (79) 

and it is an elementary exercise to show that (j6]) is the Euler-Lagrange equation for 

7r.(a,7):=sup^{n,,,(fj)}. (80) 

Van der 

WaalsESl 

interpreted the existence of two global maximizers of fl79p as a 
phase transition between a uniform gas and a uniform liquid phase of the hard- 
sphere fluid; however, since (I7j) also has uncountably many interface type solutions 
which maximize Ila^-y{ri), eventually the uniform solutions were interpreted as pure-, 
the interface type solutions as mixed phases describing the physical coexistence of 
locally pure phases. 

Our goal in this section is to prove the flnite volume analog of this gas f-)- liquid 
phase transition when the fluid is conflned in a macroscopic container A and in 
contact with both heat and matter reservoirs. Of course, the analogy can go only so 
far: with our neutral mechanical boundary conditions there are no spatially uniform 
solutions to dH]), so that the thermodynamic notion of a "pure phase" cannot apply 
in the strict sense of its original deflnition. Yet, empirically ^^3] (and intuitively) 
finite size distortions of the spatially uniform pure phases are limited to boundary 
layer effects near the container walls, so that in a macroscopic container which is 
connected to a matter reservoir the pure phases of the infinite volume thermodynamic 
formalism are approximately achieved in most of the container's interior by quasi- 
uniform density functions. On the other hand, interface type solutions will not 
maximize ^ [t]] , for the formation of an interface comes at the price of an "interface 
penalty" which becomes negligible only in the thermodynamic (infinite volume) limit. 
To be sure, we have not been able to verify all those details. What we have been 
able to prove is stated in our 



31 



Theorem 6.1: Let K he a convex container of macroscopic proportions, i.e. 0(A) ^ 
1 an(^ 0(A)/|A|i/3 = 0(1). LetV e L^{^^). Then for a subset of Q^^f{\\V\\^) at least 
three ordered hard- sphere fluid solutions of ^ exist, (at least) two of which are locally 

V stable. The extension of this subset o/ 0^i?(||^|| J to the open set 9^^ of (at least) 
triplicity of hard-sphere fluid solutions of ^ contains a phase transition curve along 
which (at least) two distinct hard-sphere fluid solutions maximize V^^[ri] globally. 
The transition is of first order in the sense of Ehrenfest, i.e. the partial derivatives 
of {a,'y) P^[a,'~f) are discontinuous across this grand canonical phase transition 
curve. 

Proof of Theorem 6.1: For the proof we adapt the hne of reasoning of Ref. [29] 
where a canonical phase transition is proved for V given by a class of regularizations 
of and p given by the perfect gas law. Yet many more technical estimates are 
needed for the current proof, which makes it somewhat long, and so we begin with 
its outline. 

In the first part of the proof we establish the multiplicity of solutions claimed in 
Theorem 6.1. We use Propositions 3.1 and 3.4 according to which a pointwise smallest 
hard-sphere fluid solution r/" ('^) of dHD exists when (a, 7) G 0«*^(|| J, and that 77™ (r) 
is locally V stable (see Prop. 3.1 in Ref.p, and also below). Also by Propositions 
3.1 and 3.4, a locally V stable pointwise largest hard-sphere fluid solution ?7"(r) of 
([S]) exists, but Proposition 3.1 left open the possibility that 77^(1") and ^^"(r) are 
identical. We will show that when A is a convex container domain of macroscopic 
proportions, then ?]"('') < V^i'^) a subset of pairs (0,7) e 0^'/(||^||J- This 
will be achieved by showing that for the favorable subset of (0,7) G 0^'f(||^||J the 
pressure functional V^^^lrj] evaluated with ^^"(r) is bounded above by a bound which 
is surpassed by the evaluation of V^^^lr]] with r/*^^. This implies that the locally 

V stable pointwise minimal solution is not a global maximizer, so another solution 
of ([8]) exists which is, yet it does not establish that this solution is a hard-sphere 
fluid solution. This in turn is guaranteed by imposing the "no non-fluid solutions 
condition" (|35i) of Proposition 3.3 on (0,7) G ©^'/^(H^H J, which leaves us with a 
bounded but non-empty set of favorable (0,7) values for sufficiently large A. This 
set is then extended by continuity to the multiplicity set introduced in Theorem 
6.1. 

In the second part of the proof we establish the existence of the phase transition 
in 6^j. Having already established, in part one, that the locally V stable pointwise 
minimal solution is not a global maximizer of (rj) when "a and 7 are big enough," 
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we recall our uniqueness results to establish that the pointwise minimal solution is in 
fact the unique global maximizer of V^^{ri) when "a and 7 are small enough." The 
rest of the proof consists in using continuity arguments to show that for favorable 
(a, 7) the pointwise minimal solution is a global maximizer of V^^lrj) but not the 
only one. 

This ends the outline of our strategy of proof. 

So our first task is to estimate V^.^l'i]'^] from above. Since each solution fj^^^ of 
dnj is a constant solution r 1— fj^^^ of (jTj), when restricted to A, this constant solution 
is a strict supersolution for ([8]) with the same (a, 7), and so the small solution of 
([H]) is necessarily bounded above by the small solution of (jS]), i.e. r]^{r) < ?/^w! see 
Proposition 3.4. Incidentally, we also know that ^J^^ — Vi uniformly for all small 
solutions of (ED when (0,7) G Qlf{\\V\\^). Also, -{V * l)^(r) < HV^H^ Vr G A. These 
pieces of information allow us to find the following upper estimate to V^^^lrjl'], 

KAv:] = / ^.(7 - i^V * d(r))rf'r + ^l^j^aVi\r- f |)C(r)C(f ) d'r <ff 

< P.h + « ll^ll, ^:dw)|A| + I ll'^^d'^ - f l)^r(r)r/r(f) d^'r £f. 

= (na,^(Cw) + |« ll^ll, Cdw) |A| + \a j^r^-{v)iy * C)A(r) d\ (81) 

and so 

iArv^,,[c] < n„,,(Cw) + 1«( ii^ii, + {v7iy * ^Da), ) , (82) 

where ( ■ )^ denotes average over A w.r.t. normalized Lebesgue measure. We will 
next show that 

ll^ll,^v:w+ {rAV*v7).), = O[0(A)-2/3]_ (83) 
We abbreviate dist(r, 9A) = sir), and llV^d ■ |)|| =11^^11 • We note the 

obvious pointwise estimate 

-{V*lUr)>\\Vl^^^. (84) 

Now we add zero, in the form of — , to r.h.s.l lH^ . then average the so 

rewritten ( 18^ over A w.r.t. normalized Lebesgue measure, multiply by the constant 
function r]^^, and get 

- > \\VlvZ^-{\\Vl - \\Vl^^^)^r]Zl (85) 
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The average at r.h.s. fl85l) is estimated as follows. The integrand - ^ ^ 

depends on r only through s(r) = dist(r, 9A), and when extended to all s > 0, it 
is decreasing fast to zero (at least like Cs~^) for s large (in molecular units); just 
asymptotically expand (A. 4) and (A. 5) for large R. Hence, and since A is convex, 

\Vl-\\Vl^^y\<C{V)\dA\ (86) 



where 

C(V) = ^°°(||V||.-||V(|-|)ll,,J<ifl (87) 
is independent of A. With fl86l) . fl87|) inserted into fl85l) . we obtain the estimate 

WlvZl < - {vZ^iV * fJ-w)J, + CiV)vZl\dA\/\A\ (88) 
with |<9A|/|A| = O(0(A)-i), by hypothesis. So for the l.h.s.dHS]) we arrive at 

\\v\\,vz^+{v:{v*v:x)^ < 

- {vz^iv * vzj^ - v:{v * v:\)^ + o[0(A)-i] = 

- <(& - V:)iV * ivZ^ + v:)\)^ + O[0(A)-i] . (89) 
The last displayed integral in ( l89l) we estimate thusly, 

- <(c.w - v:)iv * ivz^ + < 
-2((&-r/r)(v^*c;w)J, < 

-2((Cw-^r)(^*Cw).3\ = 
2ni,Cw((Cw-0), = 

^\\v\\,vZ^{vZ^-{v:\) ■ (90) 

To estimate {ti'^)j^, we recall that for (0,7) G 6^'f(||V^|| J the small hard-sphere fluid 
solution riZ < vZw — Vv that for such pairs (0,7) the map r] ^',(7 — {aV * 
r/)^(r)) with p',{-) = 5'^^(-) is convex when restricted to C°(A) fl Bjj^. Jensen's 
inequality then gives 

(r/D, >p'.(7-((«^*Oa)J- (91) 
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We now notice that 

((^*C)J, = (C(^*1)J, (92) 

and recall our pointwise estimate f lH^ and that H^H^ — 1|^|| ^ ^ is decreasing to zero at 
least like Cs~^ for s large in molecular units. So if 5A C A is a corridor of thickness 
O[0(A)-'^/^] next to the boundary (9A, then upon splitting A = (5A U (A \ 5K) we get 

(O, >P:(7^" + «r'(C)J, (93) 

where 

y = 7-O[0(A)-^/3] ^94) 

and 



t' = \\V\1-0[^{K)-\ (95) 
So (^7™)^ is a supersolution for the algebraic fixed point problem 

V = P,{l' + o^r%, (96) 

and this yields the lower bound 

(O, > (97) 

where T]^ is the smallest solution of ( p6|) . We also know that (r/X')^ ^ — Vv so 
by the concavity of r/ t— )■ (^) for f] <r]^ we easily find the explicit lower bound 

r/' > rf, (98) 

where r]^ solves the linear algebraic equation 

a.ivZj + g'.ivZJiv - VZJ = 7 - 0[(Z)iA)-'^'] + ai\\Vl - O[0(A)-i])r^, (99) 
which gives 

rf = VZ^'0[(Z)iA)-'/% (100) 

Hence, 

- (C), < O[0{A)-'/% (101) 
All in all, this proves f l83|) . i.e. we have shown that 

\M~'K,,[v:] < n„,,(Cw) + o[0(A)-2/3]. (102) 
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We next recall that the maximum of 'Pq^[?7] is estimated below (in particular) by 

^e^^K,,M>K,,[v:U- (103) 

We estimate V^^^lq^'^^] from below by using that our aV = A^Vyf + AyVy is of 
molecular effective range, and that A is convex and of macroscopic proportions, so 
that our earlier 5-corridor estimate tells us that we can find i} = 1— O[0(A)^^] G (0, 1) 
and ^ = 1 - O [0(A) -2/3] g (0, 1) so that -{V * l)lr) > ^9 Vr G ^A. Using also 
that -{V* l)^(r) < lll^ll^ Vr G A, we find 

K,MU= [ P.(7-(«^*Cw)ir))rfV + |Cw / / aV{\T-~r\)d\d'r 

J A J A J A 

> p.(7 + « lll^ll, ^9Cw)kA| - iaCw IIV^II, |A| 
= (?P.(7 + « ll^ll, ^Cw) - |«Cw \\Vl ) |A|. (104) 
Next, a simple telescoping yields the identity 

^p.(7 + « ii^iL ^^::w) = p.ii + « ii^iL c^w) - (1 - 0p.(7 + « wvivi^) - 

?(p.(7 + a ll^ll, Cw) - P.(7 + « \\V\\Jvtij) ,(105) 
and by the mean value theorem we have the further identity 

p.ii + « ll^ll, Cw) - p.(7 + « \\vi ^tiJ = (1 - WMc^ ni, Cw (106) 

for some 7 inbetween the two arguments at the l.h.s. (11061) . The monotonic increase 
of p', now gives 

p'M<p:ii + cy\\vivZJ. (107) 

The r.h.s. (11071) is independent of A and depends on (a, 7) as displayed plus implicitly 
through T]^^. Since 1 - -i? = O[0(A)"^] and 1 - <; = O[0(A)"2/3j^ conclude that 

^P.il + « IIV^II, ^C^w) > P.il + « II^IUvIw) - O[0{A)-'/% (108) 

and so 

lAr^maxP^^^M > p.(7 + a||r||^fj::j-|«Cwll^lL-O[0(A)-2/='] 

= ^a,MiJ-OmA)-'/']. (109) 
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Combining f ll02p with fll09p now yields the desired estimate 



\A\-' (^maxP^^^M - P„^,,[r/r] j > n.,,(Cw) - ^a,,{r]Zj - O[0{A)-'/% (110) 

But for (a, 7) in the triphcity set ©^'f (|| V^|| J, the A-independent function 

w{a, 7) := n,,,(Cw) - n„,,(r/-^) (111) 

vanishes only on the van der Waals gas & hquid coexistence curve a 1— )■ 7 = ■Jlt'^id), 
thereby dividing ©^'f(||^||J into two disjoint subsets, in one of which zj{a,'y) < 0, 
and zu{a, 7) > in the other. Since G7(a, 7) is independent of A, while O[0(A)~^/^] I 
as 0(A) — >■ 00, it follows that for each pair (a, 7) for which w^a,^) > 0, we have 

\A\-' (maxP^^^M - > ^(«,7) " O[0iA)-'/'] > (112) 

eventually, for large enough A. So for this subset of ©^'f (|| J, a locally V stable 
small hard-sphere fluid solution r/™ (^) < Vi exists, but it is not globally V stable. 

Our result f lll2p for the w{a,'y) > subset of ©^'f(||V^||J in sufficiently large 
A does not establish that the global maximizer is a hard-sphere fluid solution, or 
even that any other hard-sphere fluid solution of ([H]) exists for the "parameters" 
(a, 7) and A under consideration. Yet, since (11121) holds for any particular (0,7) 
in the zu{a,'j) > subset of ©^1*^(||V^|| J whenever A is sufficiently large, we can 
impose the additional condition fl35|) (with in place of $^) on (0,7), so that 

no solution to the extended exists which somewhere in A is not a (hard-sphere) 
fluid; see Proposition 3.3. Notice that fl5S]) is a sufficient but certainly not a necessary 
condition, yet to improve on it we would need to have better control over the solid 
branch of 7 1— ^,(7). In absence of such better control we consider 

(«,7)e©:'f(||V||Jn{u7(«,7)>0 A 7 + «ll^ll,^;L<7fJ- (113) 

This set, which is defined entirely in terms of the algebraic van der Waals theory with 
spatially uniform density functions, is non-empty, as can be verified by evaluating 
this van der Waals model. We conclude that when (a, 7) satisfies f lll3p and is 
fixed, then for large enough A a locally V stable pointwise minimal hard-sphere fiuid 
solution r7™(r) < rj; of ([8]) exists, but the global V^.^ maximizer is given by another. 
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pointwise larger hard-sphere fluid solution of (jS]) which is locally V stable, or locally 

V indifferent in exceptional cases. 

The existence of a third, unstable (under iterations and in V sense) solution 
sandwiched between the locally V stable minimal solution and the globally V sta- 
ble solution of ([H]) now follows via the mountain pass lemma thanks to the strong 
C°(A) differentiability of the functional V^^[rj\. By continuity we can extend the so 
constructed multiplicity sub-region of hard-sphere fluid solutions of ([8]) to a larger 
set 9^f, which is the open set of pairs (0,7) for which at least three ordered hard- 
sphere fluid solutions of ([8]) exist, (at least) two of which are locally V stable (or, 
exceptionally, locally V indifferent), and no non- fluid solution. 

This completes the part of our proof of Theorem 6.1 which establishes multiplicity 
of hard-sphere fluid solutions of ([8]) in a certain domain in (a, 7) space. We next prove 
that somewhere in this multiplicity region of hard-sphere fluid solutions a first-order 
phase transition occurs between a small and a large(r) hard-sphere fluid solution. 

We already know from part one of our current proof that for any (a, 7) satisfying 
(11131) . whenever A is large enough, then the global V^^ maximizer is given by a 
hard-sphere fluid solution of (IHl) which is pointwise larger than the locally V stable 
pointwise minimal hard-sphere fluid solution 777(1") of ([8]), which exists also, satisfies 
?7™(r) < r/j, but which is not a global V^^ maximizer in this (a, 7) region. Moreover, 
if we pick any (a, 7)0 satisfying f lll3p and pick a large enough A so that the globally 

V stable solution of dS]) is pointwise larger than rj'J^ir) for the chosen (0,7)0 and 
A, then by the continuity of the map (0,7) 1— )■ P,^{a,'^) and the continuity of the 
map (a, 7) V^^^[ri'^] restricted tcl§ r/™ < rj^, for A as chosen and now fixed, there 
exists a whole finite-measure (0,7) neighborhood of (0,7)0 in r.h.s. flll3p for which 
the globally V stable hard-sphere fluid solution of (|8]) is pointwise larger than ^^"(r), 
which in turn is not globally stable. Let this subset of r.h.s.f lTT^ be denoted by ^^O^f. 
It is a forteriori contained in the multiple hard-sphere fluid region 0^^. 

On the other hand, recall that according to Corollary 4.5 the hard-sphere fluid 
solution Tj^ of dH]) is unique if both of the following are true, pQg(7;)a$^ > 1 (with 
PcsItO ~ 0.047) and 7 < 7A(a) (with 7A(a) given in Definition 4.4). A unique 
hard-sphere fluid solution is necessarily the pointwise minimal solution, % = V"^^ 
and the conditions of Corollary 4.5 guarantee that ?7a G _^nBri{a,'y)^ so the solution 

^The map (0,7) 1— >■ "Pa^^i^yr] is generally not continuous without the size restriction on 77". For 
instance, think of the iS-shape sections of the solution diagram of the space-uniform van der Waals 
problem ([6]). 
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is "small" in the sense that 77™ <^i- If we supplement the conditions of Corollary 4.5 
with the condition 7 + a ||\^||^ r/^^^ < 7^, then no solution to the extended ([8]) exists 
which is somewhere in A not fluid, and then the unique hard- sphere fluid solution 
automatically is the unique maximizer of V^^ for any compliant (0,7) point. 

Now note that the condition on a in Corollary 4.5 is fulfilled for all (a, 7) 6^^, 
and so is the no-non-fluid condition 7 + a ||^||^ r/^^^ < 7fs. Hence we conclude from the 
discussion in the previous two paragraphs that along any constant-a ray which begins 
in *0^j and continues to arbitrarily negative 7 values there occurs a discontinuity in 
the map 7 ^ {^A'^(r)} from 7 to the set of global maximizers of V^^ which are all 
fluid. Indeed, along any such ray the constant-a map 7 1— > 77™ furnishes the unique 
global V^^ maximizer when 7 is negative enough, i.e. f]^{r) = ri1^{r) for 7 negative 

enough. Moreover, this map 7 1— )■ 77" (r) extends continuously differentiably^ into the 
region ^0^f, for which a hard-sphere fluid solution r]1'^{r) > f]J^{r) of (jH]) is the global 
maximizer of V^^, while ^^"(r) is not. Furthermore, by the local V stability of the 

pointwise minimal solutions^ a branch of pointwise non-minimal global maximizers 
cannot bifurcate off of this continuously differentiable branch of pointwise minimal 
small solution. Hence, some discontinuous change in the set of global maximizers 
must happen along each such ray, as claimed. We next clarify the nature of the 
discontinuity. 

For fixed suitable A and a as just described, we now define 7gf(a) to be the 
supremum of 7 values for which ri'^'^{r) = vTi''^) < Vi ^■^ the unique global maximizer of 
V^^^ for all-f< -f^,{a)- clearly, -^^{a) < -i%{a) < 7.t(a We also define ,7^^,(a) 

to be the infimum of 7 values for which //^ l'^) < Vi ^■^ ''^ot a global maximizer ofV^^ for 
all 7 G (*7g';(a) , + e) for some e > 0; clearly, 7^^,(0) <*7g^.(«) < 7.t(a ll'^^IIJ- 

We next show that *7gf(a) =1^1(0:). 

Indeed, suppose *7^^(a) 7^ l^eict)- Then 7^f(Q;) <=K7^^(a), and now it follows from 
the definitions of 7^^(a) and //^^(a) that ^^"(r) < fj; is a global maximizer of V^„^ 
for all 7^(,(a) < 7 <*7g^(a), though not the unique one. But then, not only are the 
values of V^^[r]] the same for its pointwise minimal maximizer 77™ and for its other 
maximizer(s) r^^, also the 7-derivatives of V^^lrj] must be the same for r^^ and for 
Tjf. By the implicit function theorem, the derivative of 7 H- V^^^lr]^] exists along any 
constant-a section of a solution branch of ([HD, except at the bifurcation points where 
it might or might not exist, but in any event either the left or right derivative w.r.t. 
7 exists, then. Now, since any currently contemplated (0,7) is not a bifurcation 
point for the pointwise minimal solution, the partial 7 derivative of V^^lrj"^] exists 
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at (a, 7). Any other maximizer of V^^[ri\ belongs to a different solution branch, and 
we may assume that in general it is not at a bifurcation point either, so the partial 
7 derivative of 7^0,7 ['7a] generally exists at the contemplated (0,7), too. Now, with 
the help of ([8]) one can easily show that, away from bifurcation points, 

= [ - * V.)i'^))d'r; (114) 
Ja 

note that in terms of the functional for the total number of particles f[T^ we can 
re-express this derivative as d-yV^^^[r]j^] = So we conclude that the two 

maximizers 77^ and r^™ of V^^It]] also have the same A^, i.e. J\f^[ri'^] = N'^[ti'J^], which 
is impossible because r^^ is the pointwise minimal solution for the given (0,7). This 
proves that *7^^(a) = 7^<,(«)- Incidentally, the proof also shows that r/™ ('^) < Vi ™^ 
a global maximizer ofV^,^ for all 7 > 7^^ (a). 

Now, by the continuity of the maps (0,7) t-)- PA{a,-y) and (0,7) 1— )■ V^^^lrj"^] 
at (a,7^^(a)) it follows that a (-)■ 7^^(a) is continuous in it's (restricted) domain of 
definition0 Moreover by the continuity of the maps (0,7) h-)- Pf^(a,'y) and (0,7) h-> 
"^0,7 ['7r] &t {(^:'lge{(^))j we also conclude that the pointwise minimal solution 77™ 
is certainly a global maximizer of V^^^ also at (a,7g^(a)), and then denoted r^^^. 
However, the definition of 7g^(a) leaves it open whether or not rj^ is the unique 
maximizer also at (a,7g^(a)), in which case the "sup" in the definition of 7g^(a) 
could be replaced by "max." We next show that at (a, 7^^(a)) the global maximizers 
of V^^^ is not unique. 

Let (0,7) be a point on a ray emanating from '^0^^ which is to the right of but 
near the curve 7 = 7^^ (a). Then, by the just proven fact that 7^^ (a) =*7^^(«), and 
by their definitions, it follows that some hard-sphere fluid solution rij^{r) > r]™ (^) of 
dH]) is a global V^.^ maximizer for each 7 > 7g^(Q;) in a right e-neighborhood of 7^^ (a). 
For each such 7 > 'y^),{c() pick such a maximizer 7]j^{r){> rj'^lr)) and consider the map 
7 (-)■ ri^{r). Notice that the set {7 > 7^^ (a)} is open. Since 7 1— )■ ^',(7) is monotonic 
increasing, each solution rj^ of ([H]) is a supersolution for ([S]) with 7 replaced by 7 — e. 
This implies that the branch 7 H- ri/y{r) of the globally stable hard-sphere fluid 
maximizers of which are bigger than rjj^ir), is pointwise monotonic increasing in 
7. We conclude that the following limit exists pointwise and strongly in C^{A)nBjj< , 

lim VA=:vt^- (115) 

^Once again, the restriction ?7™(i") < Vi is vital; without it, the domain of definition of a 7^^ (a) 
can be extended to all a € K+ but then this map is not continuous. 
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The strong C° continuity of 77 p',{'y — {ciV*ri)/^) implies that 77^^ also solves By 
the continuity of 7 Pa (a, 7) it follows that 77°^ is also a global maximizer of V^ ,y. 
This proves that at (a, 7^^(a)) the global maximizers of V^,^ is not unique, interpreted 
as a first order phase transition in the sense of mathematical coexistence of two 
distinct global maximizers of V^,^, not to be confused with the physical coexistence 
of two locally pure phases, separated by an interface, described by a single solution. 

It remains to prove that the phase transition is of first order also in the sense 
of Ehrenfest. In Ref. [33] we showed that for fixed A, the map (0,7) 1— t- Pa(«,7), 
defined in ( ITOj) . is the limit of a family of functions which are convex in a and 7, and 
so itself (bi-)convex in (a, 7), thus continuous and almost everywhere differentiable 
in both variables. We now show that at (a, 7^^(a)) generally there is a kink in both, 
a t-)- P^(q;,7) and 7 h-> P^(a,7); we ignore the exceptional case 7g^(a) = c. In that 
case by the bi-convexity of (a, 7) (-)■ P^(a,7) both partial derivatives jump up when 
crossing the grand canonical phase transition curve from smaller to larger a or 7 
values. 

By repeating almost verbatim the arguments used to prove that *7^^(a) = 7^f(a), 
one proves that the 7-derivative of Pa(q;,7) jumps at 7'\(a). Now, since we ignore 
when 7gf(«) = constant, locally there exists 7 ag^(7), the inverse function of 
7g^(a), and arguing almost verbatim again, but now using also that < 0, for 7 
suitably fixed we find that 

hm vr = vZ^ (116) 
too. Moreover, away from bifurcation points of a solution branch (0,7) ^ of 

daK,M = -k I lv{\v-i\)r^,{r)r^,{i)dhd^f. (117) 

JK J A 

Incidentally, we can also re-express this derivative in terms of the functionals for 
the total number of particles ( !T3|) and total energy (fT^ of a density function r], viz. 
c^daV^^^lr]^] = ^^f^[l]^] — S^[ri^], but we will not use this rewriting. Now, since 
(a, 7'\(a)) is not a bifurcation point for the pointwise minimal solution, also the par- 
tial a-derivative of V^^^l'i]'^] exists at (a, 7^^ (a)). Since at (a, 7^^ (a)) a pointwise larger 
global maximizer of V^^lrj] exists, too, which belongs to a solution branch which con- 
tinues to carry the global maximizers for some right neighborhood of (a, 7* (a), it 
now follows from (11171) together with V < that also the a-derivative of VAia,^) 
jumps. The transition is therefore of first order in the sense of Ehrenfest. ■ 
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The maximizer 77°^ is of pointwise minimal type fj^ir) and called the gas solution. 
The pointwise larger maximizer rjl'^ we call the liquid solution (|8]), although our proof 

does not establish that 77^*^ is quasi-uniform up to a boundary layer; numerically^^ 
this is the case, though. 

Theorem 6.1 and its proof do not establish that r^^^ is the pointwise largest hard- 
sphere fluid solution rj^'; more generally it does not establish that rj^' is the global 
maximizer of V^^lr]] for (a, 7) satisfying f lllSp and A sufficiently large. The proof 
only shows that under these conditions the global maximizer of V^^[ri] is a hard- 
sphere fluid solution which is larger than the pointwise minimal solution 77™. While 
this necessarily implies that 77" (r) > 77™(r) V r G A it does not even imply that 
"^0,7 [^A^] — 'Pa,'y[VA]- This result in turn holds in strict form and in more generality, 
as we show next. 

Proposition 6.2: // (0,7) lies in the ^^7(0,7) > subset 0/ (||V^|| J and A is 
sufficiently large so that Ul^) holds, then V^^^[ri"] > V^^^li]"^]. Thus ri^'{r) > ?7™(r) 

V r G A, and an unstable third hard-sphere fluid solution is sandwiched inbetween. 

Proof: Pick any (a, 7) in the w{a,-y) > subset of 0^'/^(||^|| J, and suppose A is 
big enough so that f lll2p holds. Consider the iteration ?7("-+^) = ^',(7— (aV^ * r/^"-*)^), 
starting from 77^°) = 77^^^ < 7/<. By Proposition 3.1 and 3.2 it iterates downward 
to the pointwise maximal solution rjff in the truncated cone C°^_(A) fl -B^jm^. As 
remarked after Proposition 3.1, a priori the minimal and maximal solutions may 
coincide, but having proved flll2p . this cannot happen because the functional V^^^lq] 
increases along any monotone converging iteration sequence. Indeed, setting 77'^"^ = 
^(n+i) _ ^(n) ^y^q difference of any two subsequent iterates, and [t/]'-"^ = (^r]^'^~^^^ + 
r/^"'')/2 for their arithmetical mean, if r]^'^\r) 7^ V r G A, then by the mean value 
theorem (apphed to the p integral) and a binomial identity (applied to the V double 
integral) we have 

P^,,[77("+^)]-P^,,[7/(")]= / {p:{j^-\T))-l4-\r))i-aV*fj^\{r)d'r > 0, (118) 

Ja 

where 7^"^(r) =7 — (aV * {[?7]}*-"'')y^(r) and {[77]} '■"''(r) is a (bounded) continuous func- 
tion which takes values between the smaller and the larger one of the two iterates 
r]^"'\r) and r]^"'~^^\r)] the inequality in f lllSp holds because —V > and, for mono- 
tonic iterations, p',(7''"'*(i")) — M'""'Hr) has the same overall sign as 77^")(r). So, 
'Pa,'y[VA] > "^0,7 [^a]' which implies that 7/*^ ^ 77™, and therefore 77" (r) > ?7™(r) 

V r G A. Thus we have established the existence of at least two distinct hard-sphere 
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fluid solutions in the w{a, 7) > subset of 6^f*^(|| ) when A is big enough so that 
ffml) holds. 

Moreover, in this limits of respective monotone iterations associated 

with the Gateaux derivative of V^^ both the pointwise minimal solution ri'^{r) and 
the pointwise maximal solution ?7"(^) locally V stable, or at most V indifferent 
in exceptional cases. Save such exceptional cases the existence of a third, unstable 
(under iterations and in V sense) solution sandwiched between ^^"(r) and rj^^ij:) 
now follows, via the mountain pass lemma, from the local V stability of these two 
solutions and the strong C°(A) differentiability of the functional V^^^l^]]- ■ 

Remark: Note that the assumptions in Proposition 6.2 do not imply that the global 
maximizer of V^^ is a hard-sphere fluid solution; for this we need to assume more, 
e.g. (I113p as in the proof of Theorem 6.1. For the proof of Proposition 6.2 we 
therefore could not assume that 77" ^ 77™, but instead had to (and did) prove it 
anew. □ 

Remark: Variants of Theorem 6.1 and Proposition 6.2 and their proofs hold with 
0« (ll^lli) replaced by G^'f ($^) and r]^^ by r/™ < as well as r]^'^^ by rjf > r/,. 
In that case the L^(M^) integrability of V is not required so that we can even add 
A^V^ to aV . If = Vn, then (jH]) can have many more than three hard-sphere fluid 
solutions for the same (a, 7), cf. Ref.[5Tl |32], even though the algebraic fixed point 
problem fl26|) has at most three solutions, still. This shows that naive inferences from 
the algebraic fixed point problem f l26|l onto the integral equation ([8]) are not to be 
drawn. □ 

Remark: Neither the proof of Theorem 6.1 nor the one of Proposition 6.2 imply 
that the pointwise maximal hard-sphere fluid solution of ([8]) is the global maximizer 
of V^^ when the pointwise minimal solution is not. If one could show, possibly 
by the index theorems of Ref. [H HO], that in the no non-fluid solutions regime at 
most two locally stable hard-sphere fluid solutions exist, but otherwise arbitrarily 
many unstable hard-sphere fluid solutions, then the pointwise maximal solution is 
the global maximizer whenever the pointwise minimal solution is not, and vice versa. 

□ 

In the (a, 7) region where the pointwise minimal solutions of ([8]) are not globally V 
stable, by their local V stability they are V metastable. Such metastability regions 
usually terminate at a spinodal line, the location of which in (0,7) space can be 
estimated. Namely, on the one hand we already know that a metastable rjj^ < rj^ 
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exists for each (a, 7) G B^j'^($]r3) fl {w{a,^) > 0} whenever A is sufficiently large. 
On the other hand, by our nonexistence result of any solution which would take only 
hard-sphere fluid values, we also know that neither a nor 7 can be arbitrarily big. 
Yet, for bounded A we can get a more subtle result, valid even if A^Vn is added to 
aV. 

Proposition 6.3: Let aV = A^Vw + AyVy + Aj^V,^, with A^, Ay, and non- 
negative, and let > be the spectral radius of —{V * ■ )^ for A C bounded. 
Assume that 

av^ > min «;(r/) = g'^{%) ^ 21.20. (119) 

r?G(0,l) 

Letrj^ =r]^{aVj^) denote the smallest solution to the equation 

av^=9'M. (120) 

and set 

7(«^a) = 92{V<) - av^V<- (121) 
Let ri{r) <f]^= (7;) ^ 0.130 be a small fluid solution of Then 7 < ^{aVj^). 

Proof: For any container A C with finite Lebesgue measure |A|, each kernel 
aV = A^V^ + AyVy + Af^Vf^, with A^, Ay, and non-negative, is a Hilbert- 
Schmidt kernel (i.e. V G L^(A x A)), and so the positive definite operator —{V* ■ )^ 
is a compact operator on L^(A). By the Krein-Rutman theorem, the spectral radius 

> of — (V * ■ )^ is the largest eigenvalue of — (V^ * ■)^, its eigenspace non- 
degenerate, and the corresponding eigenfunction nonvanishing everywhere. 

Now let ^7™(r) < rj; once again be the pointwise smallest solution of ([8]). Since 
r.h.s. (IH]) is acting as a strictly convex function on the truncated cone C^^_^_nBrj^ , we can 

apply Fujita's strategy*^ as generalized by Amann.^ Let ^(r) be the eigenfunction 
of — (\^ * ■ )^ for Vj^, normalized as probability density function so that it integrates 
to 1. Let ( • ) be the averaging functional w.r.t. ^. Taking now the average of ([8]) 
with ( ■ ) and using Jensen's inequality, we find 

which cannot be satisfied if 7 > ^{av^^). ■ 
Proposition 6.3 implies the existence of a 7^ (a) < 7('^'^a) such that 

r]ry.= lim vr\a, (123) 
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which exists and solves ([H]), satisfies the following alternative: either ||'7™JIc^(a)= Vv 
and then the constant-a section of (a, 7) t— )■ r/™ may continuously extend to 7 > 
7^(a), only then with r/™(r) ^ r/j for some r G A; or ||?7™J|cO(a) < Vv then 
the constant-a section of (a, 7) (-)■ r/" is discontinuous, i.e. left and right limits 
of the map 7 77™ U at 7^ (a) disagree. When the second alternative holds, the 
metastability region for the pointwise minimal gas-type solutions terminates at the 
curve a h- >■ 7^(0), which in this case is the spinodal curve for supersaturated vapor. 

The computation of the function 7^ (a) seems generally possible only implicitly, 
through studying the family of pointwise minimal solutions rij^{r). However, its 
upper bound j{av^) can be easily computed when the spectral radius is known. 
The latter can be computed to any desired degree of precision by monotone iteration. 
Lemma 6.4: The spectral radius of the positive operator —{V* ■ )^ is given by 

lnv^= hm iln||e(")|U2(A) (124) 

n— )-oo 

where ^("'+^) = —(^V * with ^'•^^ = 1. Moreover, it is bounded by 

-{iV*l).)<v,<\\V{\ ■ 1)11 (125) 



Proof of Lemma 6.4: The identity f ll24p is just the formula for the largest Lya- 
punov exponent (= Inf^) of our linear iteration, viewed as a dynamical system. The 
lower bound in f ll25p is obtained from the Ritz type variational principle for with 
the help of the trial function ^(r) = |A|~^/^, the upper bound by applying the sharp 
Young inequalityf^ to this variational principle. ■ 

Remark: Since the evaluation of fll24p or fll25p in Lemma 6.4 may only be feasible 
numerically for general A, the following weaker estimates are of interest, too: 

\\V(\ ■ 1)11 < IIF* 111 < ||y * 111 < \\V\\ , (126) 

II VI l^ll^i(A) - II llcO(A) - II "cO(Bn) - II 111 ' ^ ^ 

with \Bji\ = |A|. The first upper bound (= $^) is elementary. The second upper 
bound (= ||V^(| ■ 1)11 ) follows from a simple radial rearrangement inequality; 

this bound is explicitly evaluated in Appendix A. The third upper bound is again 
elementary but nontrivial only if = 0. In that case, when A in the 

sense of FisherEI^, 

then both the lower bound and upper bound in (I125P converge 
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to IIV^II^ = (AwvrV4x3 + AY47r//t2)/a, which therefore is the A hmit of the 

spectral radius. □ 

We end this subsection with the observation that in the (at least) triplicity region 
of hard-sphere fluid solutions of (jS]) where the pointwise smallest solutions ?7™(r) 
are globally V stable, the locally V stable pointwise maximal solutions ^^^^(r) are 
V metastable. Clearly, the extent of this region has a lower 7 bound because of 
Corollary 4.5, and a lower a bound because of Corollary 4.3. Moreover, since we 
imposed the sufficient condition (153]) for all solutions to be fluid, we also have an 
upper bound on a given by a < 28.9 (approximately); indeed, if this bound 
is violated by a, then for no choice of 7 satisfying (with ||\^||^ in place of $^) 
is (a, 7) G ©^'fdl^llj- Recall that this bound can be improved when better control 
over the solid branch of 7 1— >■ ^,(7) becomes available. The accurate determina- 
tion of the boundary of this metastability region is generally feasible only indirectly 
through numerical solution of the problem. Numerical solutionl^S] reveals that in 
this metastability region the fluid assumes the shape of a giant liquid drop barely 
separated from the container walls by a thin layer of vapor. This metastability 
boundary is a spinodal curve which represents the smallest giant liquid drop which 
can he contained in A given (0,7). 



7 J" STABILITY AND THE VAPOR o DROP PHASE 

TRANSITION 

In this section we discuss the thermodynamic stability of our non-uniform hard- 
sphere fluid solutions in bounded containers for the thermodynamic contact condition 
"heat reservoir," i.e. what we called T stability. 

Substituting the Carnahan-Starling approximation p,f(?j) = {g^^ o g2)(j]) i^i^o the 
entropy functional (fT5l) . one can carry out the r] integration to obtain 



S^AV] = T^''[V]- [ ^(r) 
Ja 



■ , , 3 - 27]{r) 
lnr]{r) + ' 



so for a hard-sphere fluid we have, within the Carnahan-Starling approximation, 

nH = ^M-stH- (128) 

In the following, when we speak of a hard-sphere fluid density function 77 (r) as being 
globally or locally T stable, we mean a global or local minimizer of f ll28p under the 
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constraint 



M^[7]] = N. 



(129) 



Proposition 7.1: Any hard-sphere fluid density function r]{r) which is locally or 
globally T stable under the constraint M29^) is a solution ri^{r) of ^ for the same a 
but with 7 determined by the constraint U29\) . Any globally (locally) V stable solution 
of l^j is also globally (locally) T stable. 

Proof of Proposition 7.1: Since the free energy functional (11281) is strongly C^(A) 
differentiable and coercive, its local and global minimizers satisfy the Euler-Lagrange 
equation for (11281) under the constraint (I129p . When this constraint is taken into 
account in the usual manner with the help of a Lagrange multiplier 7, viz. the "null 
functional" Af^lr]] — N is multiplied by 7 and then subtracted from ^^^[77] and rj 
then varied unconditionally, a straightforward calculation gives the Euler-Lagrange 
equation 

-^ + {aV*'n)^+g,iri) = 0, (130) 

which is precisely our ([8]) with the Carnahan-Starling approximation for the hard- 
sphere fluid equation of state. So the local and global minimizers of (I128P under the 
constraint (I129p are among the solutions of ([8]), with 7 determined by (I129p . 

As for the global J-" stability, we note that the maximum P^ (a, 7) of the pressure 
functional V^^lrj] is also given by the Legendre-Fenchel transform (sending — )■ 
7) 

P^ (a, 7) = sup \iN-F^ («, AT) I , (131) 

which, upon recalling the definition of F^{a, N), can be rephrased as the variational 
principle 

P,(a,7) = snp\^Af^[v]-TM]- (132) 

Since P^(a,7) is given by the variational principle ffTOl) . which defines the globally 
V stable solutions rj'^'^ of ([8]), it follows that each rj'^'-^ also saturates the variational 
principle (I132p — for suppose to the contrary that 7A/''^ [r/^'^] — J^alvl'^] < -^aI'^'T)' 
then 'y^f'^[r|2'^] — J^alvl'^] < "^0,7 [^a*^] 5 which we show to be impossible. Indeed, after 
some straightforward manipulations of flTBl) . given by the difference of f[T^ and f[T^ . 
using only in its reverse form (11301) . and recalling one finds that 

imv.]-J'^K] = KM (133) 
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for any solution r]^ of not just those which are globally V stable. So each globally 
V stable rj^'^ also saturates the variational principle f ll3ip . But then each rj^^ also 
saturates the variational principle for global J-" stability, with = JV'^ [vl'^] ■ 

A variation on the theme of this global stability proof gives the local J-" stability 
of locally V stable solutions rj^. The proof again uses the identity fll33p . valid for 
any solution r]^ of ([H]), but replaces P^(a,7) in fll3ip by V^^^lrj^] and the variation 
in the global Legendre-Fenchel transform (I13ip by a restriction to a neighborhood 
of 7]^. m 

Remark: Incidentally, fll3ip guarantees that P^{a,'y) is convex in 7. □ 

Remark: The infimum of under the constraint = f/'^l'i]] is generally not 

given by the Legendre-Fenchel transform of V^^[ri] (sending 7 — )■ A^). Put differently, 
^a,'y[v] ^ud J^aiv] generally not convex duals of each other. As a spin-off of this, 
the reversal of the stability conclusion in Proposition 7.1 is not allowed; i.e., not all 
globally (locally) J-" stable solutions of ([8]) are also globally (locally) V stable. □ 

Remark: When V G L^(R^) we can take the infinite volume limit. The Legendre- 
Fenchel type variational principle fll32p then becomes the thermodynamic variational 
principle ^ 

7r.f(a,7)= sup {liv)^-^ - fa[v]} , (134) 



where, for each 77 G C°(M^), 

(r/) . := hill |A|-W-M, (135) 



Mv] ■■= hm lArX-M. (136) 



For spatially uniform density functions ri{r) = r], the functional (I136P for the free- 
energy density : temperature ratio of 77 takes the simple van der Waals form 

faM = - gM - I" \\y\\,v'^, (137) 

here with the local hard-sphere thermodynamics treated in the Carnahan-Starling 
approximation. Note that the infimum of ( 1136^ under the constraint (^7)^3 = f], 
denoted 

Ua,r]):= M {U[v]\ {r])^^ = f]}, (138) 

r?eC"(K^) 
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is generally not achieved by a constant function r k-)- rj, but by a piecewise constant 
Vvdwi^) ^ C'°(M'^) (PC for "petit canonical" coming in handy), and satisfies the van 
der Waals-Maxwell formula 

U{a,f]) = CH{fM}, (139) 

where CH{ ■ } denotes the convex hull. This formula for the thermodynamic free 
energy density : temperature ratio can be rigorously obtained by taking a van der 
Waals (Kac) limit with infinitely far ranging, infinitely weak pair interactions after 
the thermodynamic (infinite volume) limit^^^^^ has been taken, see Ref.[27j for 
one-dimensional, Ref.[35] for three-dimensional systems, both with Kac interations, 
and see Refs. [19], |20] for larger classes of interactions. Formula fll39p means that in 
the thermodynamic limit the van der Waals free energy density : temperature ratio 
is itself a Legendre transform, namely the Legendre transform w.r.t. 7 of the convex 
function 7 i— )• 7r,f(a,7). This example of equivalence of ensembles at the level of the 
thermodynamic functions free energy and pressure is generally false for the finite 
volume functionals, as noted in the previous remark. For certain types of V non- 
equivalence of ensembles in van der Waals-type theory occurs even in a coupled limit 
of infinite volume and infinitely far ranging, infinitely weak pair interactions. 

□ 

In our numerical invest igat ions ^3 of the canonical non-uniform van der Waals 
theory for V = V^^ and A a ball of radius 50x~^ we found J-" stable liquid drops 
surrounded by a vapor atmosphere when a \\V\\_^ ~ 31.2. Note that for this a \\V\\_^ 
value the sufficient no-non-fluid solutions condition 7+a ^ < 7ts is violated for 
the relevant 7 values used to compute hard-sphere fluid solutions, yet solid solutions 
can nevertheless be ruled out with our more refined knowledge of the solid branch. 
We remark that the droplet solutions that we found were all situated in the (a, A^)- 
region where F^{a,N) displays the "wrong" convexity which is "jumped over" by 
the grand canonical phase transition. 

Interestingly enough, our numerical studies^^ ^ revealed that the change 
from quasi-uniform vapor state to droplet state in the canonical ensemble is not 
gradual but involves another first-order phase transition which is embedded in the 
(a, iV)-region "jumped" by the grand canonical phase transition; see also Refs. [32| 
[51] . While a complete analytical proof of all the interesting details revealed by our 
numerical studies seems futile, our next theorem does assert the existence of a petit 
canonical first-order transition between a quasi-uniform vaporous and a strongly non- 
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uniform free energy minimizer with same {a,N). It generalizes our proof - from 
regularized Newtonian interactions and simpler equation of state of the perfect gas 
to the hard-sphere equation of state and interactions which include the shorter range 
van der Waals and Yukawa interactions. To state our theorem, we recall that for 
(a, 7) = (a,7gf(a)) the functional V^^^It]] has two global maximizers in the hard- 
sphere fluid regime, the gas solution ^^^'^(r) of the pointwise minimal type ?7™(r) and 
the liquid solution r]1^{r) > 77^'^ (r). 

Theorem 7.2: Let A be a convex container of macroscopic proportions, i.e. 0(A) ^ 
1 and 0(A)/|A|^/^ = 0(1), such that a ball domain B of volume \B\ = \A\/8 is a strict 
subset of A. LetV & L^{M.^) and let a \\V\\_^ G (31— e, 31+e). Then a is in the domain 
of the map a i— )■ 7 = 7'^^(a), the grand canonical gas vs. liquid phase transition curve, 
and there exists an N^i{a) G [J^^lVA^gl^-^'^iVA^]) foi^ which two distinct solutions of 
IjM) minimize J^alv] globally under the constraint N''^[rf\ = N^^{a). The transition 
between the global T minimizers is of first order in the sense of Ehrenfest, i.e. the 
partial derivatives of (a, A^) 1— j'F^(a, A^) jump at the canonical phase transition curve 
a I— )■ A^^d(tt); provided the radial symmetric decreasing rearrangements of the two T 
minimizers intersect at a single level value. 

Remark: One of the two global T minimizers is of the pointwise minimal (given a, 7) 
type ri'^ij) and represents the supersaturated vapor phase. Our proof will suggest 
that the other one is very likely of droplet type, having a high density (liquid) core 
surrounded by a low density (vapor) atmosphere, but our proof does not conclusively 
establish the existence of such a solution type for ([8]). Numerically^^ such solutions 
do exist, and they do intersect the equal- A^ vapor solution at a single level value. □ 

Remark: The global T minimizer representing a supersaturated vapor phase is V 
metastable. The global T minimizer representing a liquid drop surrounded by a 
vapor atmosphere is V unstable. □ 

Remark: A compromise between taking the infinite volume limit A — )■ (in 
the sense of, e.g.. Fisher) and to work in a strictly finite domain is to work in 
but with the restriction that all densities are periodic w.r.t. to the 3-torus = 
M^/Z^. The canonical non- uniform van der Waals theory in was studied most 
recently in Ref.[in]. For strictly finite-range pair interactions and the equation of 
state of the lattice gas model, they proved the existence of minimizers of the so- 
called Gates-Lebowitz-Penrose free-energy functional which, when restricted to a 
single fundamental cell, look like a liquid drop surrounded by a vapor atmosphere in 
a finite container A. □ 
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Proof of Theorem 7.2: For our proof we apply the strategy of Ref.[2n] where 
a canonical phase transition of the type as stated in Theorem 7.2 is proved for V 
given by a class of regularizations of Vn and p given by the perfect gas law — except 
for the Ehrenfest part concerning the a derivative of F^[a,N), for which we follow 
Ref ■ ^30j ■ We note though that the more rapid decay of V^{r) and V^(r) with r and 
the more complicated local thermodynamics require much more delicate estimates 
in the current proof. In particular, the condition of Proposition 3.3 for the absence 
of not-all-fluid solutions is too restrictive now, so that our full knowledge of the solid 
branch §3,5(7) will be brought in. With that we now begin our proof. 

First, by evaluating the algebraic van der Waals problem for our hard-sphere 
fluid one easily verifies that when a||V^|| ~ 31 and A is large, then a is in the 
domain of 7^^ (a), the grand canonical gas vs. liquid phase transition curve described 
in Theorem 6.1. Just draw a family of straight parallel lines with slope ~ 31 into 
the second figure and note that the whole fluid triplicity 7-interval for this a value 
can be covered without intersecting the solid branch. Note that the condition of 
Proposition 3.3 for the absence of not-all-fiuid solutions is violated, though. 

Moreover, since the two V maximizers along the grand canonical phase transition 
curve are pointwise ordered, ri1^^{r) < 77°^ (r) Vr G A, we conclude that J^^lrjl^J < 
■^'^[vl'e] SO that the half-open A^-interval stated in Theorem 7.2 is not empty. 

We next recall Proposition 7.1, according to which any globally V stable solution 
of dS]) is also globally J-" stable. So in particular t]^^ = rjjf^ is a global minimizer of 
J-"^ [rj] under the constraint = A/"* [^7° 9] • This global J-" minimizer r]Jf^ = rj^'^^ is of 
the pointwise minimal (given 0,7) type riJ^{r) and situated on a fixed a-section of 
the solution branch [a, N) 1— )■ riJ^{r) of quasi- uniform small solutions (< fj^), given in 
terms of the invertible parameter representation 7 i— )■ rij^{r) and 7 h-j- A^ = J\f^[riJ^] 
for each a. This representation is well-defined because for fixed a the map 7 1— )■ 
?7"(r) is pointwise increasing and (by the implicit function theorem) continuous (even 
continuously differentiable) in the half-open 7 interval (— oo,7,f(« J] containing 
7g^(a;), where 7,f(a ||^|| J is the right limit ( FM|) of the van der Waals triplicity region 
6^'j?(||l^|| J of the hard-sphere fluid (recall our Propositions 3.1 and 3.4). This 7 

interval maps into the A^ interval (0,iV(a)], where N{a) := Af'^lfjJ^] and 17" is the 
pointwise minimal solution of ([8]) for 7 = 7,f(a ||^|| J. Note that J^^[i]'^^g] < N{a), 
by Theorem 6.1. Moreover, since by Theorem 6.1 for each a in the domain of 
7"^^ the map 7 1— ?■ ?7™(r) furnishes the unique globally V stable solution rjl^ir) for 
each 7 < 7^^ (a), by Proposition 7.1 for each admissible a as stated in Theorem 
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7.2 the map N h-> 77 "(r) then furnishes a globally J-" stable solution ^^^'^(r) for each 
N G {0,m& with <^(r) ^ r7°^(r) at iV = ArA[7^°^]. 

Furthermore, by the monotonicity of 7 t— )■ = JV^ [riJ^] and the pointwise mini- 
mality of ?7™ (for given 0,7), and by its uniqueness as solution of ([8]) for 7 < jA^a) 
(see Corollary 4.5), the fixed-a section of the branch of locally stable gas solutions 
N I— i]J^{r) furnishes the unique globally J-" stable solution for each N < A/'*[r7^], 
where 77^ is ?7"(i") for (0,7) = (q;,7a(q;)). Now let N^i{a) be the supremum over 

G (0,A^(a)] for which 1— ?7™(r)(< r/j) furnishes the unique globally J-" stable 
solution for each N < N^i{a). Clearly, N^i{a) > Af'^lvl]. We also define *N^d{a) as 
the infimum over N G (0,iV(a)] for which ?7^(r)(< f]^) is not globally J-' stable 
for each N G {*N^^{a) , ^N^^{a) + e) for some e > 0. Clearly, *N^^{a) > N^i{a). We 
show: 

(a) ,N^Ma) = at;, (a); 

(c) at N^^{a) the global J-" minimizer is not unique. 
To prove claim (a) suppose that *A^^d(«) > -^id(«)- Then from the definitions 
of *A^^d(«) and N^^iict) it follows that A^ H- ?77(r)(< Vi) is globally stable for all 
A^ G (A^^d(«) , *N^d{c()), but at least one other global minimizer exists for each 
such A^ (given a). It suffices to assume that exactly one other global J-" minimizer 
77^° exists for each such A^ (given a). But then these two minimizers of J^alv] ^^ot 
only have the same J-"^ value for each such A^, also the derivative of A^ 1— J^^IVa] 
is the same for both minimizers. Now it follows right away from (I133p that along a 
constant-a section of a solution branch of (|8]) we have 

dNJ'^K] = nV.l (140) 

where r[?7^] is the 7- value for which rj^ solves ([8]). So both hypothetical global J-" 
minimizers solve ([8]) for the same (0,7) = (a,r[?7™]), but since 77^ is the pointwise 
minimal solution at (0,7) = (a,r[?7™]), it follows that Af'^[ri^'~'] > A/''^[?7™], which 
contradicts the hypothesis that both density functions are global minimizers of J-"^ [rj] 
for the same (a,A^). This proves that ^N^^^a) = N^^^a); incidentally, the same 
type of argument also proves that A^ r7™(r)(< fj^) is not globally T stable for 
N > N^^ia). 

As for (b), to prove the first inequality we recall that by Proposition 7.1 and 
Theorem 6.1 we know that is a global minimizer of J-'^ for all N G (0, A/'^fr^^'^]). 
Now suppose that beside r]™ there exists a second global minimizer of J-"^ for some 
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A''* satisfying ^f'^[r|2] < A^* < ■^'^[VA^g]- But then, by the proof of point (a), it 
follows that (— 7- ?7™(r)(< r/,) is not globally J-" stable for each > A^*, which is a 
contradiction. This proves the first inequality in (b). 

To prove the second inequality in (b) we show that for A^ = N{a) a droplet type 
density function has lower free energy than the vapor type solution fjj^ which defines 
N{a), and by continuity this will be so also for some left neighborhood of N{a). 
Since for each A^ we will only compare densities which all integrate to the given A^, 
we can ignore the TV* functionals in J-"^ and compare 

= W / o:V{\r-r\)r]{r)r]{r)d^rd^f+ / r]{r) 

J A J A J A 

(141) 

evaluated with 77™ versus its evaluation with some droplet like density of the same 
N. 

First, let A^ = N{a). Recalling the upper bound on the gas solutions riJ^{r) < r/™^ 
where the spatially uniform van der Waals solution is for the same 7 as 77™, we have 
in particular r]]^{r) < i^Zlw We apply this bound to the interaction integral, plus use 
the estimate ((V^ * 1)a)a ^ ~ ll^lli- ^^^^ apply Jensen's inequality w.r.t. uniform 
spatial average to the (negative of the) entropy integral, noting the convexity of the 
mapx^ x\nx+x{3-2x)/{l-xf, and use that (r^^)^ = A^/|A| for all A^ G (0,iV(a)]. 
This yields the lower bound on ^^[177] given by 

IM-'Atm > -\a \\V\\XJ + + (142) 

where we wrote N for N{a). Also is a function of a, and by fllOip and (?}")^ = 
iV/|A| we have that 

Cw = ^(l + O[0(A)-^/^]), (143) 

so that up to a correction of O[0(A)^^/^], we can substitute iV/|A| for rj^^; 
other way round. On the other hand, by inserting into A"^ [rj] a trial density of the 
type "liquid drop with vapor atmosphere" which integrates to iV, we get an upper 
bound on the minimum of the reduced free energy functional for A^, given a. It 
suffices to choose a spherically symmetric trial density without atmosphere, 

^.(r) = 4x.(r), (144) 



ln?7(r) 
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where i? C A is a ball whose volume is determined by setting N /\B\ = the 
pointwise largest van der Waals solution at 7(0; J. This yields the upper bound 

V 

= ^"«l'»l)B>.l|i| + l|('°|i+(?e|^)- ("5) 

Subtracting (11421) from f ll45p and using (see Appendix A.) that 

{{V*l),)^ = -\\Vl{l-O[l/0{B)]), (146) 

and anticipating that |A|/|_B| = 0[1] so that we can neglect the O[1/0{B)] correction, 
we find that the upper bound (I145p on the infimum of A is lower than the lower bound 
(11421) on the free energy of the gas solution at = when 

|A| _^ 3-2jV/|i?| _ 3-27V/IAI 

.||V||.>2 (-^l«0'^(^^ , (147) 

\B\ |A| 

up to a correction of O[0(A)^^/'^]. The criterion (I147P can be re-expressed as 



ln5s^ + 



M \^ 7" Jim 



a||y||^>2 y^H^ ^'''""^ . (148) 

vdW 

up to a correction of O[0(A)~^/^]. Now, for a = 31 as stipulated in Theorem 

7.2, the ratio 17"^ : Cdw ~ 9> with fjl'^^ ^ 0.41 and Cdw ~ 0-045. These values yield 
r.h.s.([IlS])^ 28.75 < 31, and also |A|/|S| ?a 9 > 8 so that B fits into A, satisfying the 
hypothesis of Theorem 7.2. This proves that the droplet-type density function has 
a lower free-energy : temperature ratio than the quasi-uniform vaporous solution of 
the same A^ = N{a) for a \\V\\_^ = 31. By continuity the regime where droplet type 
densities have lower free-energy : temperature ratio than the quasi-uniform solutions 
extends to an open neighborhood of the chosen a for the corresponding N[a). 

Second, by continuity again, the same conclusion also extends to some open left 
neighborhood of N{a) for each such a in the neighborhood of the chosen a || = 31. 
This completes the proof of claim (b). 



54 



Continuity and closedness arguments for the solution curves prove claim (c) in a 
similar fashion of reasoning as used in the proof of Theorem 6.1. Here we also use 
Proposition 7.1, according to which the solution rjj^ is also locally J-" stable for each 
a and N in the domain of the map t— )■ r/^ because the pointwise minimal (given 
(a, 7)) solutions rjj^ are locally V stable (we ignore the exceptional cases when 1]^^ is 
locally V indifferent). By its local J-" stability, no bifurcation off of this gas branch 
occurs for < N{a), in particular not for = A/'^fr^^^] ( < N{a)). 

Henceforth we will write r^™ = 77^^ for the quasi-uniform, vaporous J-" minimizer, 
and tJa ^ ?7a d fo'^ ^^le non-quasi-uniform minimizer of (presumed) droplet type; we 
say "presumed" for, strictly speaking we haven't shown that it is a droplet, although 
our above proof and the numerical evidence^^ suggests it is. We remark that even 
though in our proof we worked with a trial droplet without atmosphere, any droplet 
type minimizer of J-'^ [77] must solve (|8]) and therefore must have a low-density atmo- 
sphere, as r.h.s.dH]) is bounded away from 0. 

Finally, we show that the canonical vapor versus droplet transition is of first 
order in the sense of Ehrenfest, for which we need the hypothesized, yet empirically 
suggested, level intersection property. 

First, in our proof of point (a) above we showed that the constant-a derivatives 
of 1-^ J^^IVa] at A^ = A^vd(«) cannot be the same for the quasi-uniform minimizer 
rjj^l and for the non- quasi- uniform minimizer rj^'^^. So our proof of point (a) above 
already proves that the constant-a derivative of N F^{a,N) is discontinuous 
at A^ = N^i{a). In fact, the constant-a derivative of N F^{a,N) jumps down 
when A^ increases. This follows from fll40p and the monotonicity properties of the 
pointwise minimal solution branch of ([8]). For suppose that r[?7j'^] > r[?7^'!y]. Then, 
since 77^*^ ^ 77"^, which denotes the pointwise minimal solution at (a, r[77^'^]), we have 

^""[vl^] > -^HvZd] > -^""[vl^], which contradicts the fact that M^^lvl^] = -^^K!;]- 
Note that for this part of our proof of the Ehrenfest property we did not need to 
invoke that the two global J-" minimizers intersect only at a single level value. 

Next, to prove that the constant- A^ map a J^oIVaI ^ kink at a^^^^N), where 
N a = a^i{N) is the local inverse function to the curve a N = N^^{a), 
and which exists locally unless the latter is constant, we adapt, and improve on, 
the strategy of Ref . |3U] . First, being J-" minimizers, the free-energy : temperature 
ratio is the same for the quasi-uniform minimizer rj^l and for the non-quasi-uniform 
minimizer 7]^^^, i.e. J^a[v^%] = -^alvl^^]- Recalling the identity ( I133P we see that this 



55 



implies the equality 

N{r[vZ] - nvZ]) = ^:,r[.pc] [vZ] - Kn,ic, KfJ, (149) 

where r[?7^] is the 7- value for which r]^ is a solution of Now suppose that, for 
given N, the derivative of a t— )■ J^^iv^] at a = a^^{N) is the same for both minimiz- 
ers. By the implicit function theorem these derivatives exist along the constant- 
sections of the solution branches of ([8]); at a = a^^^N) the derivative for 77^^ may 
have to be read as right-derivative, though generically it will be a derivative. Using 
the variational principle for F^{a,N) and the definition of J^^lv] obtaiiil^ 

daJ'M = \ 11 - i\)v,ij)r,,ij) d'rd'f. (150) 

J A J A 

Incidentally, for our V < f ll50p shows that a F^{a,N) is monotonic decreas- 
ing, but we won't need that. By f ll50p we conclude that the hypothesized equality 
^aJ^aivl^v] — ^aJ^aiv^d] a niiuimizers at the same (a, A^) implies that 

their potential energy : temperature ratios are the same, too. Inspection of the defi- 
nition ^ of the pressure : temperature ratio functional of a density function 77 now 
reveals 

- = / b«(C.(r)) - vM,%^Wr, (151) 

J A 

where 'p,M) — 9iiv)- By inserting fll5ip into fll49p we obtain the equahty 

JA 

Recall that at the end of the first part of the Ehrenfest proof, i.e. of the disconti- 
nuity of the constant-a derivative of J^aVft^]i showed that the l.h.s. (11521) < 0. 
We now complete our proof that the constant- derivative of a J^^Wi^] is dis- 
continuous at the canonical phase transition curve provided the radial symmetric 

®This formula may cause some temporary consternation, for a thermodynamic free-energy : tem- 
perature ratio should satisfy the "thermodynamic relation" daF^[rjf^] = S^VIa] along the globally 
F stable solution branch of ([S]), with £^ [77^] given in (HH). This puzzle is resolved by noticing that 
in our strictly classical setup we have omitted even the minimal amount of quantum mechanics 
normally injected into classical statistical mechanics with the help of the de Broglie wavelength, as 
per "normalization" of the entropy and chemical potential; see our Appendix C. 
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decreasing rearrangements of the two T minimizers intersect at only a single density 
value X by showing that r.h.s. (11521) > under this provision. 

For any density function //(r) G C^(A), let *//(|r|) denote its radially symmetric 
decreasing equimeasurable rearrangement supported in the ball B of volume |-B| = 
|A|. Then, using that 

/ /(r/(r))c?V = / /(>(|r|))dV (153) 

JA JB 

for any continuous function / : R — )■ R, and using the mean-value theorem (with 
sandwiched between the two J-' minimizers uniquely determined), and invoking 
the hypothesized level intersection property of the two minimizers, we find 

Ja Jb 

= [ p:M\r\)r&\)-vZ(\^\)]d'r 

Jb 

> P'M [ rvZil^l) vZi\^\)]d'r 
Jb 

= 0, (154) 

where the inequality is readily proved by estimating the penultimate integral sepa- 
rately on the positive and negative parts of its integrand. This already concludes 
the second part of the proof of the Ehrenfest property, but we supplement the result 
by showing that the constant- derivative of a i— )■ F^{a,N) jumps down when a 
increases. 

Indeed, since for a < a^i{N) the vaporous solution is the unique global J-" mini- 
mizer, suppose now that at a = a^^i^N) we have daJ^^[ri^'i^] < daJ^^lrj^^]- -^^^ then 
by straightforward adaptation of our proof of the discontinuity of daF^{a, N) we 
now conclude that 

> iv(r[<:=j - rivZ]) = P^^^j^pej [vZ] - ViZ] > o, (iss) 

and so the derivative d^F^^a, N) must jump down at a = a^^i^N). 

The proof of Theorem 7.2 is complete. ■ 

We end section VII with some comments regarding the proof of Theorem 7.2. 
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Remark: Since our proof of the Ehrenfest property relies on the provision of the 
level intersection property of the two global T minimizers, it seems prudent to have 
a backup strategy just in case the provision turns out not to hold for non-spherical 
containers; also for spherical containers it hasn't been proven yet, although in that 
case there is numerical evidence in its favor. [33] The following argument does not 
rely on the provision of Theorem 7.2, and could be completed with some sharper 
estimates. 

Namely, we use that r^^^ < rj, is the pointwise minimal solution for (a, 7) = 
(a, r[?7^^]), so that for the same (a, 7) we have the bound 77^^ < ^()- More- 

over, rf^^ is quasi-uniform in the sense that it is nearly constant except for a small 
boundary layer near 9A, viz. (recalling ( llOip ) 

r/-w-(^r), <O[0(A)-^/% (156) 

where {t]^'^y)^ = N/\A\ is the uniform mean over A. Furthermore, since both J-" 
minimizers have equal "mass" A^, we have the identity 

(Oa = (157) 

Since x (-> p,f(a;) is a positive, increasing, convex function, dividing l.h.s.f ll5^ by |A| 
and applying Jensen's inequality combined with these identities and estimates yields 

(p.,(C.(r))>^-(p.,(C.W)>A > P.MvZ\) - P.^ivZj 

> -O[0(A)-2/3], (158) 

where the small error, due to the boundary layer effects, goes to zero as A goes to 
M^, but is not identically zero. 

Thus, to complete this proof one would need to show that the difference r[?7^*^] — 
r[r7^^] < stays away from zero; alternatively, the proof would be completed if one 
could control the error term in Jensen's inequality to the effect that (p,f(^7^*rf(i')))A ~ 
P»f ((^Id)A) > C* > independently of sufficiently large A. □ 

Remark: In the limit of vanishing hard-sphere volume the local thermodynamics 
goes over into that of the perfect gas. In this case the hard-sphere pressure : tem- 
perature ratio as function of f] is simply the identity map, and then r.h.s. (11521) is 
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identically zero, and the proof of the discontinuity of the constant-A^ derivative of 
a ^ ^aVff!^] is complete, then. This in fact is the proof of Ref . [30] . □ 

Remark: We note that the jumping down of the constant- derivative of a t— )■ 
F^{a,N) at a = a^^iN) also implies (for our V < 0) that £^[r]l%\ > ^aivl^]^ which 
is seen by noting fllSOp and recalling the definition f|T^ of the energy, keeping in 
mind the constancy of A^'^[?7^*^] at a = a^i{N). With the jumping down of the 
energy : temperature ratio at a = a^^^N), the constancy of ^a[?7pc^] at a = a^i{N) 
then in turn implies that i5^f [?7^^] > '5^f['7^^], i-e. the entropy jumps down also. □ 

Remark: Our proof of the canonical phase transition reveals two metastability 
regions in its (a, A^) neighborhood in which locally J-" stable solutions of ([8]) exist. 
Also these metastability regions should terminate at their spinodal lines. We have 
to leave the determination of their location in [a, N) space for some future work. 

□ 
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APPENDIX 



A. THE INTERACTION INTEGRALS IN SPHERICAL GEOMETRY 

In spherical geometry we are in the position to obtain several explicit results. 
Lemma A.l: Let A = Bfi be a ball of radius R centered at the origin. Then 
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Setting r = in (A.l), (A. 2), and (A. 3) produces 
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Integrating (A.l), (A. 2), and (A. 3) over yields 
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(A.6) 
(160) 

, (A.7) 

(A.8) 

(A.9) 
(161) 



B. ASSOCIATED PARTIAL DIFFERENTIAL EQUATIONS 

For Yukawa and Newton kernels the integral equation (|H]) for the density ri{r) 
in A C is equivalent to a semilinear ellipticES PDE of second order together 
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with consistent boundary condition for the corresponding chemical self potential per 
particle, —V * rj. Thus, setting — Vy * r] = ip, we find from ([8]) that ip solves 

- A^(r) = 47rp'.(7 + aij{r)) - K^tp{v) .(B.l) (162) 

In the formal Newtonian limit /t — )■ we have Vy — ?■ Kj, and (B.l) reduces to 

- A^(r) = 47rp',(7 + a^p{r)) .(B.2) (163) 

In the low density limit, (B.l) reduces to 

- Atp{r) = 47rCe"^('') - ^^^(r) , (B.3) (164) 

and (B.2) to 

- A^(r) = 47rCe°'^(") , (B.4) (165) 

with ( = e'^ the fugacity. In each case, (E]) evaluated at dA provides a nonlinear 
and nonlocal boundary condition for ip, which makes it quite difficult to study these 
PDEs in general domains. 

Remark: If in one replaces ^',(7) by the strictly convex function exp(7), then 
the alternative stated after Proposition 6.3 ceases to exist and the map 7 1— )■ ?7™|a G 
C°(A) actually terminates at 7^(a); see Refs. [T8| |5] . □ 

For spherically symmetric solutions, i.e. ip{r) = (f){r) in a ball of radius i?, 
satisfying the regularity condition 0'(O) = 0, the PDEs (B.l), (B.2), (B.3), and (B.4) 
simplify to ODEs with nonlinear and nonlocal boundary conditions that read, for 
(B.l): 

g-K/J r-R 

(f)(R)=47T—- rsmh(Kr)p',h + a(f)(r))drJB.5) (166) 
kR Jo 



for (B.2): 



for (B.3): 



for (B.4): 



1 

cf){R)=47T- p',{-r + a(f){r))dr {B.6) (167) 

R Jo 



-kR pR 

(j){R) = ^7tC—^ / r sinh(Kr)e"'^('')rfr , (B.7) (168) 
kR Jq 

1 

(j){R) = ATTC- rV'^('')cir.(B.8) (169) 

R Jo 
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For spherical symmetry, (B.2) has nice scahng properties which facihtate its dis- 
cussion and aid in its numerical integration on a machine; see Ref.[5Tl [32]. Its 
low density limit (B.4) becomes the homologously invariant isothermal gaseous ball 
equation, which has been extensively studied by Emden,^^ Chandrasekhar^I^ and 
others. 1^ ^26] g^c^ scaling properties are not shared by (B.l), or (B.3), for which 
numerical studies of radial solutions apparently have not yet been carried out. 

The spherical version of (jS]), with V = ]Ajv in does not seem to reduce to 
an ODE, and numerical integration of the integral equation (jH]) with V = are 
correspondingly more involved, see Ref . [33] . 

The special case A = M^, with rj{r) solving ([7]), is of interest in itself, as explained 
in the introduction. Since equations (B.l), (B.2), (B.3), and (B.4) do not depend 
on A, the same PDEs cover the case A = M^. However, instead of taking the formal 
limits A — M'^ for their self-consistent boundary conditions, the situation is more 
subtle. We illustrate this with the spherically symmetric situation, with finite R 
boundary conditions (B.5), (B.6), (B.7), and (B.8). In fact we need to drop (B.6) 
and (B.8), for their limits are infinite because the respective equations (B.2) and 
(B.4) do not possess solutions with their right-hand side in L^(]R^). For (B.l) under 
spherical symmetry, we may or may not include the limit -R — i- oo of (B.5), which 
is easily shown to be zero because p' is bounded; similarly, for bounded radially 
symmetric solutions of (B.3) the limit i? — ?■ oo of (B.7) vanishes. Yet if we do 
include the condition that ip{r) — )■ as |r| — )■ cxd, then we throw out all the constant 
solutions r i— )■ ijj^^^{r) = —fj^^^ IlK-Hj- This shows that the spatially constant van 
der Waals densities r i— )■ fj^^^ are more subtle limits of the finite volume non-uniform 
van der Waals densities, namely in the sense of supnorm convergence on the members 
of any sequence of nested compact subsets of M'^, which sequence converges to M^; of 
course, convergence is also weak, i.e. pointwise. 

In the wide interface approximation,^^^^^^ for our short ranged Vw G L^(R'^) 
and Vy € L^(R^) the convolution V * t] given by for A = can be expanded to 
second order, and the fixed point equation ([7]) reduces to a PDE for rj (not ip), viz. 

- aM2{V)A7]{r) + a \\V\\^r]{r) = p'.-\vir)) - 7, (B.9) (170) 

where 

M2iV) = - [ |x|V(|x|)c/V(B.10) (171) 
6 Jr3 

is the "second moment" of V. Notice that (B.9) is the Euler-Lagrange equation for 
a so-called Cahn-Hilliard functional, studied recently in Ref. [5] . 
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Our numerical studies of (|H]) for A = 5/j with i? = 50 and x = 1 revealed that 
near the critical point one finds solutions with inhomogeneity scale R. This leads us 
to the following (mildly vague) conjecture: 

Conjecture: For (a, 7) in some droplet neighborhood of the (weakly A-dependent) 
critical point, the wide interface approximation becomes asymptotically exact, in the 
sense that one finds droplet solutions rj^ ^ of ^ which converge in a suitable but 
reasonable sense to solutions of (Q.^) in some "universal" limit as A — )■ . 
C. REVERSAL TO THE DIMENSIONAL QUANTITIES OF PHYSICS 

As conventional in chemical physics we have used dimensionless units in which 
the "density" rj is actually the volume fraction occupied by all the microscopic balls. 
Thus, if N balls, having volume |6| each, are inside a container A of volume |A|, then 
7^ = A^|6|/|A|. Also, we have absorbed several "constants of nature" in our quantities, 
and moreover ignored the usual heuristic injection of quantum mechanics as per the 
thermal de Broglie wavelength. To make contact with physics one needs to reconvert 
our dimensionless into dimensional variables. It suffices to do the conversion for the 
model with the van der Waals interaction potential; the conversion for the model 
with Yukawa or Newton interactions is done entirely analogously. 

Thus, |6| is dimensional (a volume), and we have to make the following replace- 
ments "dimensionless"—)- "dimensional" quantities: r — t- r/|6|^/^ for the position 
vectors, and therefore all lengths — in particular, x — )■ |6|^/^x; next, a — > 13a for 
the coupling constant : temperature ratio; 7 — )■ — ln(A^3/|6|) for the chemical 
potential per particle : temperature ratio; p — )• \b\l3p for the pressure : temperature 
ratio; r] — )■ \b\p for the particle density. We now have {3 = {kBT)~^, with T the 
temperature in degree Kelvin, and Ajb is the thermal de Broglie wavelength. In the 
same vein, we need to replace ln?7(r) — )■ In (^(r)/^^^) in the entropy functional, where 
p^g = /h^ is the thermal "de Broglie density." 

For applications to, say, fiuids made of the nobel elements, the physical ordering 
is \b\ < 47rx"^/3 and x0(A) ^ 1. Numerically, \h\ iA^, and x~^ 2A seem rea- 
sonable, while 0(A) ^ 10 — lO^cm seems a reasonable range of laboratory container 
sizes. Also, the dimensional van der Waals coupling constant a has physical dimen- 
sion of energy, numerically in the range of "typical molecular binding energies" of the 
natural gases, although of course there is no quantum mechanical formation of Ne2, 
Ar2, etc. molecules in nature. The attraction between Ne, Ar, etc. atoms is man- 
ifested most dramatically through the condensation / evaporation phase transition 
exhibited by these chemical elements of matter. 



63 



D. ERRATA FOR REFS.pS] and 

All corrections to our previous papers Refs. |29] and [33] are easy to make. Those 
in the categories typo and slip of pen are just listed without comment, those which 
deserve a commentary are commented on in footnotes. Expressions to be replaced 
are surrounded by quotation marks. 

D.a. Errata for Ref.fE^ 

p. 223: above (2.66), replace "it has been shown*^^)" b}|§ "it has been argued^^-*". 
p.238: in (4.7aL replace "|?7" by "|f/|". 
p. 248: replacS 

"Then also s{pi) = s{p2)- That implies that there exists an incompress- 
ible mapping pi i— )■ p2- (Note that entropy is conserved for incompressible 
mappings.) As has been shown in ref. 10 (see also refs. 8 and 21), any 
given po can be mapped incompressibly to a unique spherical minimizer 
Pm of e(p) with s(po) = s{pm)- By construction both pi and p2 minimize 
e(p) under conservation of entropy; hence, pi = p2, in contradiction to 
the assumption that the densities are not identical." 

"Since also /(pi) = /(P2) (where /(p) is given in (3.28), here with ip = 0), 
we then conclude that J^exp(— /3trt^ * pi)cPr = exp(— /3trf^ * p2)d^^' as 
well. Hence, 77(^1; Pi) = ''?(Ar; P2) = ^tr (see p. 238 for the definition of 
?7(/3; p)). But then, since pi = pp^^, so that -/3trf/ *Pi = -f^tvU *p^,^ = ^^t^ 
is the unique pointwise minimal solution of (4.6) for this value of 77 = rjtr, 
it follows that exp(— /3trf^ * pi)d?r < exp(— /3trf^ * P2)d^r, which is a 
contradiction." 

^In fact, there is a small mistake in ref. 1 of Ref. |29) to the effect that the factor {N — 1) in 
(2.66) (quoted from ref. 1) is incorrect. The correct factor is A''; see Ref.|31|. 

^•^The critical sentence is: "That implies that there exists an incompressible mapping pi 1— >■ p2-" 
While true for some types of phase transitions associated with symmetry breaking, it is not clear 
that such incompressible mappings exist in the context of the theorem. M.K. is grateful to Elliott 
Lieb for kindly pointing this out. 

-'^-'^ Notice that the correction given here not only avoids the pitfall of the original proof, it also 
eliminates the requirement of the original proof that A be spherical. This nonspherical argument, 
taken from Ref. [30] . is a special case of the argument in our proof of Theorem 7.2; see the penultimate 
remark in section VII. 
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D.b. Errata for Ref.JS^ 

p. 1353: replace "The Hilbert space" by "The space" 
p.1364: in (6.19), replace "O[ro/i?]" by 'V[{ro/RfY 
p.1376: in Ref.33, replace "Phnys." by "Phys." 
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